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ABSTRACT 

An  experiment  to  measure  the  muon  capture  rate  in  liquid 
hydrogen  was  carried  out  using  scintillation  counter  techniques. 
The  experiment  tests  muon-electron  universality  in  interaction 
with  a  proton  at  a  momentum  transfer  of  100  MeV/c.  The  use 
of  ultra-pure  liquid  hydrogen,  a  purified  muon  beam,  neutron- 
gamma  ray  discriminating  detectors  and  oscilloscope  photography 
made  possible  the  detection  of  the  5.2  MeV  neutron  from  the 
relatively  improbable  capture  reaction  ( |x  +  p  -•  n  +  v)  . 
Observation  of  the  well  known  pion  capture  reaction 
(tt-  +  p  -•  n  +  <y)  served  as  a  check  on  many  aspects  of  the 
experimental  arrangement  as  well  as  on  the  neutron  detector 
efficiency  calculation.  Interpretation  of  the  muon  capture 
experiment  depends  on  knowledge  of  the  muon  molecular  wave 
function  and  the  muon-proton  spin  configuration;  these  are 
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known  from  detailed  calculations.  The  muon  molecular  forma¬ 
tion  rate  was  measured  in  an  independent  experiment,  but  the 
present  result  is  relatively  insensitive  to  it.  The  stability 
of  the  muon  molecular  ortho  state  against  transitions  to  the 
para  ground  state  for  times  comparable  to  the  free  muon  life¬ 
time  is  confirmed.  The  experimental  result  of  the  capture 
rate  measurement,  464  ±  42  sec  ^  is  compared  with  562  ±  60  sec  } 
the  rate  expected  on  the  basis  of  the  conserved  vector  current 
theory,  muon-electron  universality  and  the  accepted  value  of 
the  "induced  pseudo-scalar"  coupling  constant. 
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I .  INTRODUCTION 

A.  Status  of  the  Universal  Fermi  Interaction 

The  concept  of  a  Universal  Fermi  Interaction  (UFI) 

governing  all  weak  processes  was  placed  on  a  quantitative 

14 

foundation  by  the  discovery  that  the  rate  of  0  beta  decay 
could  be  used  to  predict  the  muon  decay  rate  rather  accurately . ^ 
It  is  of  great  interest  to  test  this  hypothesis  in  the  third 
class  of  strangeness -conserving  weak  interactions,  the  muon- 
proton  interaction. 

The  muon-proton  interaction  at  low  energies  is  available 

for  study  in  nuclear  muon  capture,  analogous  to  K-electron 

capture.  This  capture  can  take  place  because  negative  muons 

stopping  in  matter  will  eventually  cascade  to  the  lowest 

2 

Bohr  orbit  of  a  muonic  atom,  a  structure  that  is  more  tightly 

bound  than  a  normal  atom  because  of  the  large  muon  mass 

(m  =  207  m  ) .  The  atomic  wave  function  overlap  is  suffi- 
M-  e 

ciently  large  to  make  the  muon  capture  interaction  competi¬ 
tive  with  the  free  muon  decay.  Most  muon  capture  experiments 
carried  out  with  complex  nuclei  have  as  a  prerequisite  detailed 
knowledge  of  the  nuclear  wave  functions,  if  accurate  infer¬ 
ences  about  the  muon-proton  coupling  constants  are  to  be  made. 

A  long  series  of  muon  capture  experiments  in  complex 
nuclei  has  provided  some  qualitative  information  on  the 
capture  reaction.  The  axial  vector  coupling  (7)  in  muon 
capture  was  shown  to  be  equal  to  that  in  beta  decay  to  within 

4 

20-30%  by  studying  the  (purely  Gamow-Teller)  reaction  cycle 
+  C12  -  B12  (J  =  1  ground  state)  +  V,  B12  -  C12  +  e"  +  V. 
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5a  19 

A  spin  dependent  capture  rate  has  been  observed  in  F 

Measurements  of  the  angular  distribution  of  neutrons  from 

5b  5c  .  16 

muon  capture  m  calcium  and  of  muon  capture  in  0  popu¬ 
lating  discrete  levels  of  N^6  yield  results  that  are  consis¬ 
tent  with  a  large  contribution  from  the  so-called  induced 
pseudo-scalar  coupling  (see  Section  II) .  The  muon  capture 

3 

rate  in  He  has  been  observed  to  be  within  about  10%  of  that 

5ci 

expected  on  the  basis  of  the  current  theory. 

Experiments  were  initiated  to  study  muon  capture  in 
hydrogen,  a  system  free  of  nuclear  complications,  to  provide 
a  test  of  the  universality  of  the  V-A  weak  interaction  coup¬ 
ling  theory.  A  complete  experimental  result  would  indicate 
whether  or  not  a  muon  behaves  exactly  as  an  electron  does  in 
interactions  with  a  proton.  Beta  decay  experiments  involve 
momentum  transfers  much  lower  than  the  100  Mev/c  available 
when  a  muon  is  captured;  an  unknown  momentum  dependence  of 
the  effective  coupling  constants  could  result  in  an  ambigurty 
of  interpretation. 

Previous  experiments  on  muon  capture  in  hydrogen  have 
established  the  relative  sign  of  the  vector  and  axial  vector 
coupling6  and  have  established  the  presence  of  a  strong 

7 

vector  contribution  . 

A  muon-proton  coupling  scheme  of  the  form  vector  minus 
axial  vector  (V-A)  results  in  a  strong  spin  dependence  of  the 

g 

muon  capture  reaction 

\l~  +  p  -*  n  +  v  (1) 

If  the  reaction  proceeds  from  a  singlet  state  then  the  neutrons 
are  emitted  isotropically  and  are  in  fast  polarized  due  to 
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the  two  component  nature  of  the  neutrino.  On  the  other  hand 

3 

if  the  muon  and  proton  are  in  an  L  =  0  triplet  state  (  S) 
then  the  reaction  is  inhibited  altogether.  Momentum  depen¬ 
dent  contributions  to  the  coupling  modify  the  V,A  destructive 
interference  and  allow  some  triplet  capture  to  take  place, 
the  singlet  rate  being  altered  as  well. 

B.  Experimental  Problems  in  Muon  Capture  in  Hydrogen 

1.  Specification  of  the  Initial  State 

The  possibility  of  a  spin  dependent  reaction  rate  makes 

accurate  specification  of  the  initial  muon  proton  state  vital 

to  an  interpretation  of  experimental  results.  The  history 

2 

of  a  negative  muon  after  considerable  reduction  of  velocity 
in  liquid  hydrogen  (density  =  0.07  gm/cc)  proceeds  as  follows 
(see  Fig.  1) .  At  a  low  velocity  the  muon  is  captured  by  a 
proton  and  cascades  to  th°  lowest  Bohr  orbit  (the  Is  state) 
in  about  10  sec,  this  atomic  system  is  formed  in  a  statis¬ 
tical  mixture  of  triplet  and  singlet  states.  The  (y,p)  system 

is  small  (Bohr  radius  =  a  =  a  •  m  /m  ) ,  neutral  and  com- 

|i  jb  e'  p," 

pletely  detached  from  the  hydrogen  molecule  from  which  it 
was  formed.  It  is  able  to  wander  through  the  hydrogen,  the 
muon  being  exchanged  from  proton  to  proton.  If  the  liquid 
hydrogen  is  pure  the  result  of  these  exchanges  is  a  100% 
population  of  the  energetically  lower  singlet  state  and  total 
depolarization  of  the  muons .  If  even  small  amounts  of  deu¬ 
terium  or  other  atoms  are  present  then  it  is  more  or  less 
likely  that  one  of  the  exchanges  will  leave  the  muon  bound 
to  a  deuteron  (or  heavier  atom)  which  it  will  not  subsequently 
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leave.  In  the  case  of  (|j,d)  the  reduced  mass  effect  makes 
that  system  more  tightly  bound  by  130  eV  than  a  (pp)  system. 

In  pure  liquid  hydrogen  the  singlet  (p.p)  system  will  eventu¬ 
ally  (in  about  0.6  (isec)  become  part  of  a  (p^ip)+  molecular 
ion^  if  the  muon  has  not  decayed  first  (with  lifetime  2.2  |j.sec)  . 

The  molecular  ion  is  a  stable  configuration  and  it  is 

from  this  system  that  the  muon  is  most  likely  to  be  captured 

or  to  decay.  In  this  experiment  the  time  of  each  capture 

event  is  measured;  thus  if  the  transition  rate  X  from  a 

PP 

(M-P)  system  to  a  (ppp)  is  known  then  the  configuration  from 
which  capture  takes  place  is  also  known.  The  transition  rate 
has  been  measured  in  a  separate  experiment .10,11  Qniy  capture 
events  which  have  originated  in  a  (p|j,p)  system  have  been 
included  in  the  analysis.  Measurements  of  the  muon  capture 
rates  using  a  bubble  chamber  include  contributions  from  both 
atomic  and  molecular  capture  since  time  measurement  is  impos¬ 
sible.  There  are  two  questions  concerning  the  nature  of  the 
molecular  ion  which  must  be  considered  and  which  have  been 
answered  by  detailed  calculations;  (1)  What  is  the  angular 
momentum  character  of  the  molecular  ion  and  hence  what  is 
the  composition  in  terms  of  singlet  and  triplet  (p,p)  ? 

(2)  What  is  the  wave  function  overlap  of  the  muon  and  proton? 

The  (p(ip)+  molecular  ion  is  almost  invariably  formed  in  the 
ortho  state  (para  formation  is  <  10  smaller)  and  it  has 
been  shown  with  a  high  degree  of  confidence  that  only  states 
with  total  spin  =  1/2  are  present.  It  then  follows  (see 
Appendix  IV)  that  the  capture  rate  from  the  molecular  ion 
is  given  by  2y  (3/4  X  +  1/4  X  ),  where  the  Vs  refer  to 

OS  t 
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the  capture  rates  in  the  triplet  and  singlet  atoms  and  2yq 

is  a  wave  function  overlap  factor  for  the  ortho  molecule. 

Although  the  para  state  is  the  molecular  ion  ground  state, 

it  can  be  shown  that  transitions  to  it  are  extremely  unlikely 

12 

during  the  muon  lifetime .  The  capture  rate  m  this  state 
is  lower,  however,  being  given  by  2y  (1/4  X  +  3/4  X ), 

P  S  u 

and  any  appreciable  amount  of  conversion  to  it  could  be 
detected  in  this  experiment. 

2 .  The  Expected  Rate 

On  the  basis  of  the  V-A  interaction  with  certain  cor¬ 
rection  terms  included,  the  rate  of  reaction  (1)  when  the 

g 

muon  is  bound  in  the  molecular  ion  is  expected  to  be  about 
560  sec  ■*".  The  decay  of  the  muon 

|jf  -  e"  +  V  +  V  (2) 

which  is  not  appreciably  affected  by  the  muon  being  bound 
in  a  hydrogen  orbit  (binding  energy  =2.5  keV)  takes  place 
at  the  rate14  0.455  x  10^  sec  1 ,  so  that  capture  is  only 
1/1000  as  likely  as  decay.  In  the  case  of  nuclei  of  higher 
atomic  number  the  capture  rate  increases  at  least  as  fast 
as  Z  due  to  the  increased  muon-nucleus  overlap  (aQ  ~  l/z) 
and  is  equal  to  the  decay  rate  at  Z  =  12  (magnesium) . 

3 .  Detection  of  Capture  Events 

The  only  signature  of  a  muon  capture  event  is  the  detec¬ 
tion  of  the  5.2  MeV  neutron  in  reaction  (1).  The  experi¬ 
mental  arrangement  for  measuring  the  muon  capture  rate  must 
be  able  to  distinguish  this  neutron  from  all  background  and 


must  detect  it  with  known  efficiency  relative  to  the  number 
of  available  muons.  Because  of  the  possibility  of  transfer 
of  a  muon  to  deuterium  or  other  impurities  in  the  hydrogen 
the  experiment  must  be  carried  out  under  ultra-pure  conditions. 
Whereas  in  a  gas  the  various  molecular  transfer  rates  would 
be  under  control  of  the  pressure  and  it  would  be  possible 

p 

to  alter  at  will  the  ratio  of  triplet  and  sirglet  states 
beam  intensity  considerations  make  it  difficult  to  use  a  gas 
target  at  reasonable  pressure,  thus  liquid  hydrogen  was 
chosen  as  the  target  material.  Gamma  rays  produced  in  the 
bremsstrahlung  of  decay  electrons  [reaction  (2)]  constituted 
an  intense  source  of  neutral  background  which  could  not  be 
eliminated  by  shielding.  The  neutron  detectors  that  were 
employed  were  capable  of  distinguishing  between  neutrons  and 
gamma  rays  eliminating  much  of  the  background  problem;  they 
were  simulaneously  used  for  decay  electron  detection. 

The  known  rate  of  the  pion  capture  reaction 

TT+p-n  +  y  (3) 

emitting  an  8.9  MeV  neutron  made  it  possible  to  perform  a 
control  experiment  in  situ. 

The  experiment  to  be  described  here  is  an  extension  and 

refinement  of  a  similar  experiment  previously  performed  at 
7 

this  laboratory.  Evaluation  of  the  earlier  experiment 
resulted  in  minor  improvements  in  the  target,  meson  beam, 
and  running  procedures  that  were  incorporated  into  the  present 
experiment. 
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II.  THEORETICAL  BACKGROUND 


The  muon  capture  rate  is  given  by 


w 


1  ksE 
+ 


(4) 


where  k  is  the  neutrino  momentum,  a  is  the  muon  atomic  Bohr 
radius  and  M  is  the  matrix  element.  This  has  been  derived 
by  assuming  a  hydrogenic  wave  function  for  the  muon  initial 
state 

Y(r)  =  (Tra3)"^  exp  (-r/a)  (5) 


and  then  integrating  the  matrix  element  over  all  momenta  that 
are  relevant;  the  integral  can  be  approximated  by 

M  =  (2tt)3/2  M(0)  f(0)  (6) 


where  M(Q)  is  the  matrix  element  evaluated  at  zero  initial 
momentum  and  Y(0)  is  the  wave  function  evaluated  at  the  origin. 
The  matrix  element  can  be  written^5  as 

M  ~  2  ^i^v°iU|j.  <njJi|P>  (7) 


where  the  Ch  are  the  usual  operators  S,  V,  A,  T,  P  and  the 
sum  is  taken  over  as  many  of  them  as  are  present.  In  the 
limit  of  small  momentum  transfer  and  treating  nucleons  as 
though  strong  interactions  were  absent  the  matrix  element 
reduces  to  the  form 


M  ~  l  9i(uv°iV  <  V>iV  +  gi("v°iuji)  ("n°iY5up)  +  h'c-  (8) 

If  in  addition  the  interaction  is  taken  as  V,A  and  the 
coefficients  of  the  odd  parity  terms  g|  in  the  Hamiltonian 
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are  equal  in  magnitude  to  those  of  the  even  parity  terms, 
the  matrix  element  becomes 


"~sv5vyi+vy"»vv  +gAuvV5(i+vv5nVy  +  h-°- 


(9) 


In  the  non-relativistic  limit  the  vector  operator  reduces 

—4 

to  unity  while  the  axial  vector  operator  reduces  to  a ,  the 
spin  matrix  so  that  the  matrix  element  becomes  proportional 
to 


«  ~  gv  +  gA  <  <y  V 


(10) 


The  term  (a  •  ap>  can  be  evaluated  from  the  relationship 
(a^+cr.p)2  =  3  +  3  +  2 o,,‘ for  the  singlet  and  triplet  cases 


X  ~ 

s 


gv  -  3g 


M-  p 
2 


and  \ , 


9v  +  9a 


(11) 


Using  the  relationship  (0^*0^)  +  2(a^-Op)-3  =  0,  the  square 

of  the  matrix  element  can  be  written  in  general  as 


M2  ~  gv2  +  3gA2  +  2(gvgA  -  g*2)^'^) 


(12) 


The  table  below  lists  numbers  (M  ;  proportional  to  the 


capture  rate  for  either  of  the  two  initial  spin  states; 
we  assume  that  the  coupling  constants  are  equal  in  magnitude. 
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TABLE  I 


Muon  Capture 

Matrix  Element 

Initial 

State 

Coupling 

Const. 

( a  •  a  ) 

M.  P 

M2 

singlet 

gV  =  gA 

-3 

4g  2 
yV 

triplet 

II 

+1 

<? 

to 

singlet 

gV  =  -gA 

-3 

!6gv2 

triplet 

gV  =  -gA 

+1 

0 

Thus  the  equality  of  sign  of  the  coupling  constants  results 

in  no  hyperfine  effect  while  the  expected  V-A  theory  results 

in  a  much  larger  capture  rate  in  the  singlet  state  but  no 

triplet  state  capture  at  all. 

In  order  that  the  Hamiltonian  be  invariant  under  time 

reversal  the  coupling  constants  must  be  relatively  real. 

In  the  case  of  the  free  neutron  decay,  the  phase  angle  between 

g^  and  g^  was  shown"^  to  be  approximately  180°.  If  it  is 

assumed  that  | g^  |  =  j  g^  |  but  that  they  have  a  relative  phase  cp, 

then  the  expressions  for  the  singlet  and  triplet  muon  capture 

rates  become  TTg  ~  gv  (10  -  6  cos  cp)  and  X^  ~  gv  (2  +  2  cos  cp) 

and  a  system  involving  3/4  singlet  capture  and  1/4  triplet 

2 

would  then  have  a  capture  rate  X  ~  g^  (8-4  cos  cp)  .  The 
usual  choice  of  cp  =  v  results  in  the  maximum  capture  rate. 

The  above  treatment  of  the  strongly  interacting  particles 


is,  of  course,  not  legitimate.  The  part  of  the  Hamiltonian 
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which  concerns  the  neutron  and  proton  can  be  handled  phenomo- 
logically  by  writing  the  vector  and  axial  vector  part  of  the 
matrix  element  (n|j\ |p)  in  the  most  general  form  consistent 
with  the  known  conservation  laws.  Assuming  Lorentz  invar¬ 
iance,  parity  invariance  and  time-reversal  invariance,  which 

.  17 

hold  for  strong  interactions ,  we  write 

<n|jv|p>  =  %  {Vv(q2)  +  CT\piSFM(q2)  +  iq\FS(q2)}  UP  •  (13) 

(n|JAjp)  =  Gn  {YxFA(q2)  +  aXpiqpFpT(q2)  +  iqxFp(q2)|  7^.  (14) 

The  J  quantities  are  currents  with  the  appropriate  isotopic 

spin  properties  to  convert  the  initial  proton  state  to  the 

final  neutron  state,  the  u's  are  Dirac  wave  functions;  the 

entire  complication  of  strong  interactions  has  been  absorbed 

into  the  F's  which  are  form  factors  evaluated  at  the  momentum 

transfer  of  the  reaction  (q  =  0.95  m  ) .  At  zero  momentum 

|l 

transfer  the  above  expressions  reduce  to  the  familiar  form 
with  Fv(0) ,  Fa(0)  equal  to  1  and  the  other  terms  vanishing 
because  q  =  0.  The  functions  of  momentum  transfer,  F,  can 
be  interpreted  as  amplitudes  for  interaction  diagrams  involving 
intermediate  pion  states.  See  Fig.  2  for  examples  of  such 
diagrams . 

The  requirement  that  each  term  in  the  matrix  elements 

iTTT2 

have  the  same  properties  under  G  transformation  (G  =  Ce  ) 
allows  two  terms  to  be  eliminated:  Fg  and  Fpip.  The  hypothesis 
of  the  Conserved  Vector  Current  holds  that  the  vector  part 
of  the  interaction  Jv  is  a  conserved  quantity  and  that  it 

V 

can  be  obtained  from  the  vector  electromagnetic  currect  Jx 
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by  a  rotation  in  isotopic  spin  space the  weak  interaction 
current  being  proportional  to  T^  -  iT2  and  the  other  to  T^- 
As  a  consequence  of  this  assumption  one  can  identify  the  two 
terms  FT7  and  Fw  as  follows: 


Fv(q2)  =  Fx(q2) 

(15) 

FM(q2)  =  (l/2M)(pp-nn)  F2(q2), 

(16) 

where  F1  and  F2  are  the  electric  and  magnetic  form  factors 

for  nucleons  and  the  jx's  are  anomalous  magnetic  moments: 

2  2  2 

jXp  =  1.793  and  jj_n  =  -1.913.  Using  also  F^(qQ)  for  qQ  =  0.9  m^ 
the  matrix  element  becomes  proportional  to  the  following 
expression  (v  is  the  neutrino  momentum) 

M  ~  Gv  +  ga  <VV  •  GP  <(V  V)(vv))  (17) 

where  the  new  coupling  constants  are  defined  as 
GV  =  gV^1  +  2M^ 

GA  =  gA  gV^X  +  ^p  -(Xn*  2M 

GP  =[gP  ~  gA  ‘  V1  +  ^p"  ^n}]2M  (18) 

In  this  expression  the  coupling  constants  gv<  gA  and  gp  are 
proportional  to  F^ ,  FA  and  Fp;  they  are  the  muon  dressed- 
nucleon  coupling  constants  associated  with  muon  capture. 

The  assumption  of  a  universal  Fermi  interaction  in  which 
muons  and  electrons  are  coupled  in  the  same  way  in  weak  inter¬ 
actions  implies  that  the  "bare"  coupling  constants  that  are 
effective  in  beta  decay  and  muon  capture  are  equal.  The  muon 
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dressed  coupling  constants  and  the  electron  dressed  coupling 
constants  then  differ  because  of  momentum  transfer  differ¬ 
ences.  The  relationship^  between  these  two  sets  of  coupling 
constants  is: 


g£  =  0.97  g£  (19) 

g£  =  0.999  gP  (20) 


The  induced  pseudoscalar  coupling  constant  gp  has  been  calcu¬ 
lated,  using  dispersion  theory  and  the  pion  decay  rate,  and 

18 

has  been  found  to  be  equal  to 

9p  =  8gA  .  ( 21) 

19 

Beta  decay  experiments  on  the  free  neutron  and  the  decay 

14  8  8 

of  0  permitted  evaluation  of  the  ratio  as  soon  as  it 

was  clear  that  the  Fermi  and  Gamow-Teller  matrix  elements 

contain  no  contribution  other  than  V  and  A.  Experiments  give 

g£  =  -1.2  g$  (22) 


The  muon  capture  rate  in  hydrogen  can  finally  be  written 

Q 

explicitly  if  the  measured  beta-decay  rate  is  used  to  eval¬ 
uate  the  numerical  constants. 

-1 


uu 


=  158  [(gj)  2  +  3(gJ)  2]  (a  +  b<yap>)  (23) 


where 


2  2 

3  -  V  +  3rV 

b  -  2<Sga  -  ?  GVGP>  -  2<ga2  -  !  gagp> 
rA2  *  GA2  +  J  <GP2  -  2GAGP> 


(24) 
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Wit  h  the  accepted  values  of  the  coupling  constants,  one 
has  for  the  singlet  and  triplet  atomic  capture  rates: 

X  =  636  sec  ^ 
s 

=  13  sec 


A  calculation  including  higher  relativistic  corrections  and 

a  slightly  different  value  for  the  ratio  of  beta  decay  coup- 

8  8 

ling  constants  (g^/g^  =  -1.25)  yields  somewhat  different 

.  .  20 

capture  rates : 


X  =  713  sec 
s 


-1 


X  =  16.9  sec 


-1 


The  calculated  values  of  the  muon-proton  overlap  factors 

3  12 

(in  units  of  l/fra^  )  for  the  ortho  and  para  molecules  are 

27  =1.165  27  =1.308. 

o  p 

Due  to  uncertainties  in  the  molecular  wave  function,  these 

21 

numbers  are  estimated  to  be  correct  to  within  about  10% 

(the  wave  function  expansion  coefficient  is  m^/M^) .  The 
muon-proton  spin  correlation  in  the  molecular  states  (see 
Appendix  IV)  implies  molecular  capture  rates: 

gA/gV  =  -1'2  gA/gV  =  ~1-25 

Xpara  =  27p(^Xg  +  f^)  =  220  sec  1  249  sec  1 

560  sec  ■*"  630  sec  ^ 


Xortho  =  27o(fXg  +  iXt) 
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There  have  been  three  previous  measurements  of  the  hydrogen 
capture  rate.  One  was  performed  in  this  laboratory  and  was 
the  precursor  of  the  experiment  herein  described;  the  others 
were  both  performed  using  a  bubble  chamber  to  detect  the  stop¬ 
ping  muon  and  the  neutron  induced  proton  recoil.  In  the  bubble 
chamber  experiments  capture  events  are  recorded  which  took 
place  in  both  the  atomic  and  molecular  systems  and  the  expected 
capture  rate  (using  Xpp,p  =  1.9  x  10  sec  )  is  then  about  4% 
higher  than  that  in  the  molecule.  The  three  previous  results 
are 

7  -1 

Columbia  515  ±  85  sec 

Chicago  22  435  ±  100  sec-"*- 

CERN  23  420  ±  75  sec-1  . 

III.  EXPERIMENTAL  PREREQUISITES 

A.  Target 

Because  of  the  possibility  that  a  muon,  even  after  having 
reached  the  ground  state  of  a  p,p  atom,  may  be  transferred  to  an 
impurity  of  higher  atomic  number,  it  is  essential  to  carry  out 
the  muon  capture  experiment  in  ultra-pure  hydrogen.  Assuming 
the  transfer  rate  to  impurities,  such  as  nitrogen,  to  be 
5  x  10"^  sec-^  (several  times  that  to  deuterium)  and  a  concen¬ 
tration  of  one  part  per  million  (1  ppm)  then  the  number  of  muons 
ending  up  bound  to  the  impurity  atom  is  2^%  of  the  number 
available.  If  the  muon  capture  rate  were  100  times  greater  in 
the  impurity  than  in  hydrogen  the  effect  that  one  seeks  to 
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measure  would  be  overwhelmed  by  nitrogen  captures.  It  is  then 
clear  that  one  must  take  extraordinary  precautions  with  the 
purity  of  the  hydrogen;  even  the  presence  of  deuterium  cannot 
be  tolerated  at  concentrations  greater  than  a  few  parts  per 
million.  Isotopically  pure  hydrogen  was  purchased  from  l'Aire 
Liquide,  Paris.  This  hydrogen  was  further  purified  at  the 
entrance  to  the  hydrogen  target  by  passing  it  into  a  heated 
palladium  coil;  only  hydrogen  forms  the  appropriate  chemical 
bond  with  palladium  that  allows  it  to  diffuse  through  the  walls 
of  the  tube ,  other  gases  do  not  penetrate  the  coil  and  are 
eliminated.  This  purification  system  was  subjected  to  a  series 
of  stringent  tests  of  which  two  examples  follow.  In  a  test  for 
nitrogen  transmission  the  system  following  the  palladium  purifier 
was  evacuated;  nitrogen  under  pressure  was  admitted  to  the  input 
end  of  the  palladium  coil  and  the  output  pressure  was  watched. 

No  change  in  pressure  was  detected,  implying  a  maximum  leak 

3 

through  the  palladium  of  1  part  nitrogen  to  10  parts  hydrogen. 

To  test  for  helium  transmission  a  conventional  helium  leak  detec¬ 
tor  was  connected  to  the  output  end  of  the  purification  system 
and  the  input  end  was  placed  under  helium  pressure;  no  net  in¬ 
crease  in  helium  at  the  output  was  detectable.  The  closed 

“6 

system  "leaked  up"  at  the  rate  of  5  x  10  mm  of  Hg/hr.  at  a 
temperature  of  50°  C. 

Liquid  hydrogen  was  chosen  as  the  target  material  for  the 
muon  capture  experiment;  its  specific  gravity,  while  only  0.07 
(near  atmospheric  pressure) ,  is  much  greater  than  that  of  the 
gas  at  convenient  pressures.  The  low  temperature  of  the  liquid 
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(20.4°  K)  should  retard  the  outgassing  of  foreign  materials  from 
the  walls  of  its  container.  Because  the  experiment  must  be 
carried  out  in  ultra  pure  hydrogen,  the  gas  phase  was  purified 
and  then  liquefied  within  the  target  itself;  this  was  accom¬ 
plished  by  keeping  the  target  region  closed  and  in  thermal 
contact  with  a  reservoir  of  normal  liquid  hydrogen.  At  a  very 
slight  pressure  above  atmospheric  the  entering  gas  from  the 
purification  system  liquefied  in  the  previously  evacuated  target 
volume  (see  Fig.  3) .  The  filling  procedure  took  approximately 
20  hours  during  which  time  the  normal  liquid  hydrogen,  absorb¬ 
ing  the  heat  of  condensation,  was  frequently  replenished. 

Once  the  target  was  full  heat  losses  were  quite  small.  The 
target  remained  full  for  several  weeks  with  no  deterioration 
in  purity  and  required  no  maintenance  except  for  refilling  of 
the  liquid  hydrogen  and  liquid  nitrogen  reservoirs. 

Before  filling,  the  "pure  system"  was  evacuated  to  a 
pressure  of  10  mm  and  baked  for  several  days  to  release 
materials  adsorbed  on  the  walls.  An  ion  pump  backed  by  a 
chemical  adsorption  pump  was  used  to  evacuate  the  system; 
this  is  preferable  to  the  more  usual  oil  diffusion  pump  since 
it  eliminates  the  possibility  that  carbon  and  carbon  compounds 
will  contaminate  the  volume  being  evacuated.  At  liquid  nitro- 

-7 

gen  temperature  the  pressure  was  6  x  10  mm  of  Hg. 

The  necessity  of  subjecting  the  target  to  a  wide  range  of 
temper at ure  and  pressure  conditions  dictated  to  some  extent 
the  choice  of  target  materials.  The  primary  requirement  was, 
however,  that  it  be  constructed  exclusively  of  elements  with 
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high  atomic  number  to  minimize  the  effect  of  muons  which  might 

stop  in  the  target  walls .  Although  the  muon  capture  rate 

3 

increases  roughly  as  Z  ,  the  muon  lifetime  in  the  atomic  ground 
state  is  already  as  low  as  0.2  p.sec  in  the  case  of  iron  as  com¬ 
pared  to  2.2  (isec  in  hydrogen  so  that  one  has  only  to  delay  data 
taking  for  a  )isec  after  which  the  entire  effect  of  wall  captures 
has  disappeared.  The  inner  hydrogen  target  was  constructed  of 
stainless  steel  for  mechanical  strength  and  was  lined  with  sil¬ 
ver  to  provide  a  very  high  Z  barrier  for  muons  headed  toward 
the  walls.  A  silver  thermal  radiation  shield  held  at  liquid 
nitrogen  temperature  was  suspended  inside  the  vacuum  jacket 
whose  outer  wall  was  fabricated  of  zirconium.  Figure  4  shows 
the  thickness  and  location  of  the  target  walls  and  window.  The 
target  window  through  which  the  muon  beam  enters  was  made  as 
thin  as  safety  permitted  in  order  to  decrease  the  number  of 
muons  stopping  in  it  and,  more  important,  to  minimize  multiple 
scattering  of  the  beam  which  would  increase  the  likelihood  of 
wall  captures  in  the  vicinity  of  the  neutron  detectors . 

Table  II  lists  the  composition  of  the  target  side  walls 
and  front  window.  Under  the  conditions  of  the  experiment  with 
60  MeV  incident  muons  and  26  MeV  muons  entering  the  target,  one 
can  estimate  the  root  mean  square  scattering  angle  due  to  mul¬ 
tiple  scattering  in  the  front  window  and  in  the  2f  in.  of 
beryllium  moderator  which  precedes  it.  The  front  window  produces 
an  rms  angle  of  4.5°  while  the  moderator  scatters  through  roughly 
2.3°.  Since  the  multiple  scattering  in  the  hydrogen  is  quite 
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negligible,  at  least  until  the  very  end  of  the  muon  range,  one 
can  predict  the  rms  spread  of  the  beam  in  the  vicinity  of  the 
neutron  detectors  due  to  these  two  sources.  The  beryllium  is 
at  a  mean  distance  of  16  in.  from  the  neutron  detectors,  whereas 
the  front  target  window  is  only  9  in.  from  the  detector  center 
line.  Thus  the  spread  due  to  scattering  is  0.6  in.  due  to  the 
moderator  and  0.7  in.  due  to  the  window. 

It  is  difficult  to  use  these  numbers  to  estimate  the  number 
of  wall  stopping  muons ,  partly  because  the  beam  entering  the 
moderator  may  itself  have  an  angular  divergence  of  perhaps  8°, 
considering  the  collimator  geometry. 

The  pure  hydrogen  chamber  in  the  target  extends  for  a 
distance  of  13  in.  beyond  the  neutron  detector  center;  the  back 
wall  is  zirconium  but  it  is  extremely  unlikely  that  any  muon 
can  penetrate  the  entire  length  of  the  target.  The  "dirty" 
liquid  hydrogen  cooling  reservoir  and  the  liquid  nitrogen  supply 
which  cools  the  radiation  shield  are  situated  behind  the  pure 
target  which  protrudes  forward  into  the  detector  assembly. 

B.  Meson  Beam 

While  the  hydrogen  target  was  in  preparation,  studies  were 
undertaken  to  improve  the  composition  of  the  available  meson 
beam.  In  normal  operation  of  the  Nevis  Synchro-cyclotron  the 
circulating  proton  beam  strikes  a  target  located  at  a  radius  of 
about  72  in.  Negative  pions  produced  in  the  forward  direction 
are  bent  sufficiently  in  the  fringe  field  of  the  cyclotron  to 
pass  through  an  evacuated  channel  in  the  main  shielding  wall 
(see  Fig.  5) .  Momentum  selection  is  accomplished  by  the  action 
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of  the  fringe  field  and  the  small  acceptance  angle  of  the  beam 
channel.  The  only  muons  which  are  able  to  enter  the  experimental 
target  area  are  those  which  originate  from  pion  decays  in  the 
vicinity  of  the  cyclotron  target;  the  typical  radius  of  the 
muon  source  is  2  in.  or  3  in.  Most  muons  from  decays  along  the 
flight  path  strike  collimator  walls  and  are  lost.  The  negative 
meson  beam  is  also  rich  in  electrons  which  are  produced  when 
the  gamma  rays  from  the  decay  of  neutral  pions  materialize  in 
the  cyclotron  target.  The  usual  particle  ratio  in  the  negative 
beam  is  approximately 

Pions  :  7 

Electrons  :  7 

Muons  :  1  . 

For  the  purpose  of  performing  the  muon  capture  experiment  this 
situation  was  considered  unsatisfactory  for  the  following  rea¬ 
sons:  First,  pions  stopping  in  the  vicinity  of  the  hydrogen 

target  would  produce  uncorrelated  neutron  background  with  a 
wide  energy  spectrum;  Second,  electrons  in  the  beam  would  be 
electronically  indistinguishable  from  muons  and  would  make  an 
accurate  determination  of  the  number  of  muons  stopping  in  the 
hydrogen  quite  difficult;  they  would  also  increase  the  number 
of  accidental  delayed  "muon"  -neutron  coincidences. 

Since  the  source  of  muons  is  more  diffuse  than  that  of 
pions  and  electrons  it  is  possible,  by  suitable  adjustment  of 
the  bending  magnet  current  outside  the  shielding  wall,  to 
reduce  greatly  the  transmission  of  the  unwanted  particles  while 
at  the  same  time  to  decrease  the  number  of  muons  in  the  beam  by 
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less  than  a  factor  of  2.  Searches  for  the  optimum  conditions 

for  this  purified  beam  were  undertaken  using  an  analog  time-of- 
28 

flight  system.  Two  scintillation  counters  (#0  and  #2)  employ¬ 
ing  56AVP  photomultipliers  were  placed  about  6  meters  apart; 
the  output  pulses  from  these  were  used  as  the  start  and  stop 
pulses  in  a  1/4  nsec  resolution  time-to-height  converter  whose 
output  was  displayed  in  a  ion  channel  pulse  analyzer.  At  the 
particle  momenta  of  interest  (about  120  MeV/c)  the  pions  and 
muons  could  be  easily  resolved.  The  mass  spectrum  that  was 
obtained  in  this  way  was  checked  from  time  to  time  by  measuring 
the  differential  stopping  ..ue  of  the  beam  as  a  function  of 
absorber  thickness.  The  w.riing  conditions  that  were  finally 
chosen  involved  an  increac  :■  in  bending  magnet  current  to  select 
muons  of  128  MeV'c,  this  r  ntum  is  about  15%  higher  than  the 
momentum  of  the  majority  of  the  pions  which  are  normally  trans¬ 
mitted  through  the  shielding  wall  (beam  channel  "D") .  Under 
these  conditions  less  than  2%  of  the  beam  particles  were  pions 
and  electrons. 

Typical  differential  beam  stopping  distributions  as  a 
function  of  moderator  thickness  are  shown  in  Fig.  6  for  the 
purified  and  unpurified  situations.  Fig.  7  shows  the  mass 
spectrum  of  the  bearr:  consisting  of  pions,  muon  and  electrons 
as  the  purification  is  improved. 

In  normal  cyclotron  operation  the  proton  beam  strikes  a 
stationary  target;  for’ about  500  jisec  the  protons  circulating 
with  increasing  radius  continue  to  produce  mesons .  Shortly 
after  the  ejection  time  the  radio-frequency  accelerating  voltage 
is  turned  off  and  the  cyclotron  is  dormant  until  the  frequency 
modulating  condenser  blades  have  rotated  into  proper  phase  for 
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the  beginning  of  the  next  acceleration  cycle.  Thus  only 

1/2  msec  out  of  the  16§  msec  period  is  used  for  experiments. 

The  cyclotron  "duty-cycle"  is  said  to  be  equal  to  30  in  this 

case,  the  average  time  between  mesons  being  30  times  shorter 

than  it  would  be  if  they  were  produced  at  a  uniform  rate  at  all 

times.  Close  spacing  increases  the  chances  of  a  random  or 

accidental  coincidence  between  a  muon  pulse  and  a  background 

neutron  pulse  also  produced  during  ejection  time.  A  method 
24 

was  developed  to  improve  the  duty-cycle  by  a  factor  of  10. 

The  fixed  target  is  replaced  by  a  target  which  vibrates  in  a 
vertical  plane  in  phase  with  cyclotron  operation.  The  accel¬ 
erating  voltage  is  removed  when  the  proton  beam  reaches  the 
radius  of  the  target;  it  continues  to  circulate  at  this  radius 
while  the  target  moves  upward  into  the  beam.  Vertical  oscilla¬ 
tions  of  the  proton  beam  together  with  the  relatively  slow 
motion  of  the  target  permit  beam  ejection  during  the  entire 
8  msec  that  the  radio-frequency  accelerating  voltage  is  off. 

The  experimental  electronic  logic  is  gated  on  only  during  this 
time  and  receives  about  75%  of  the  total  beam  available  with 
fixed  target  operation. 

C.  Background 

In  view  of  the  very  low  yield  of  neutrons  to  be  expected 
from  the  muon  capture  process ,  extensive  measures  were  taken 
to  reduce  the  neutron  background.  A  directional  neutron  detec¬ 
tor  was  constructed  by  embedding  a  cylindrical  2  liter  volume 
of  liquid  scintillator  capable  of  neutron-gamma  ray  discrimina- 


tion  in  a  5  ft  cube  assembly  of  steel  and  lead,  which  was 
slung  from  a  crane  over  the  cyclotron  experimental  floor.  The 
array  had  an  angular  resolution  of  20°  to  30°,  and  was  capable 
of  moving  horizontally  and  vertically  as  well  as  being  rotated. 
Several  weak  points  in  the  shielding  wall  were  located  and 
corrected,  and  the  plane  of  proton  orbits  was  verified  to  be 
a  source  of  neutron  background.  As  a  result  of  these  studies 
an  optimum  working  location  on  the  cyclotron  floor  was  chosen 
and  a  shielded  "house"  was  constructed  of  6  in.  thick  iron 
barbettes;  extra  steel  and  zinc  shielding  was  installed  in 
the  direction  of  the  open  channel.  It  was  decided  to  arrange 
for  the  neutron  detectors  used  in  the  experiment  to  be  located 
out  of  the  median  or  proton  plane.  Once  the  major  sources  of 
background  had  been  eliminated  the  remainder  was  observed  to  be 
rather  isotropic  in  angular  distribution.  A  steel  roof  sev¬ 
eral  feel  thick  was  placed  over  the  basic  structure  and  the 
interior  of  the  "house"  was  lined  with  paraffin.  An  effort 
was  made  to  use  diverse  materials  for  shielding  to  avoid  trans¬ 
mission  of  neutrons  at  certain  specific  energies  due  to  pecu¬ 
liarities  of  individual  neutron  cross  sections.  In  the  vicinity 
of  the  neutron  detectors  boron-loaded  paraffin,  which  lacks  the 
annoying  2.2  MeV  capture  gamma  ray  characteristic  of  most 
hydrogenous  materials,  was  used;  finally,  several  inches  of  lead 
surrounded  each  neutron  detector.  As  a  result  of  the  extensive 
shielding  arrangement  the  background  counting  rate  was  less 
than  0.25  neutrons  per  second  per  liter  of  detector  in  the 
energy  range  of  interest.  The  effective  background  rate  was 
reduced  further  by  electronic  requirements  to  be  discussed  below. 
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D.  Measurement  of  Molecular  Constants 

Using  the  same  hydrogen  target  and  beam  conditions ,  an 

experiment  was  performed  to  measure  several  molecular  constants 

as  well  as  (pd)  fusion  parameters.  Of  greatest  interest  to 

the  present  experiment  was  measurement  of  X  ,  the  ( |xp)  to 

rP 

(PM-P)  conversion  rate.  This  is  covered  in  greater  detail  in 

11 

the  thesis  of  E.  Bleser,  and  will  be  reviewed  briefly  here. 

In  the  presence  of  a  small  amount  of  deuterium  some  of 

the  (|ip)  atoms  are  diverted  from  their  normal  tendency  to  form 

the  (p^,p)  molecular  ion.  The  muon  which  is  continually  being 

exchanged  from  proton  to  proton  may  become  attached  to  a 

deuteron  from  which,  because  of  its  increased  binding  energy, 

it  will  no  longer  be  transferable  back  to  a  proton,  but  will 

form  (pp.d)  which  may  undergo  fusion.  A  measurement  of  the 

3 

yield  of  the  5.5  MeV  fusion  gamma  rays  (p  +  d  He  +  y)  as  a 

function  of  the  deuterium  concentration  (c)  determines  X 

PP 

if  the  transfer  rate  to  a  deuteron  (X  )  is  known.  One  method 

e 

of  extracting  the  number  is  as  follows:  The  total  yield  of 
fusion  gamma  rays  is 


kcX 


X  +  X 

o  pp 


+  cX 


as  can  be  seen  considering  the  three  available  channels  for  a 
(i.p  atom;  muon  decay,  formation  of  (p|i,p)  and  formation  ^d.  The 
multiplicative  constant  k  is  the  gamma  ray  yield  when  c  is  large, 
(the  saturation  value  of  the  yield) .  The  derivative  of  this 
expression  in  the  limit  of  low  concentration  is 
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X  + 
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It  turns  out  that  X  is  several  times  larger  than  X  so  that 

PP  o 

the  slope  is  in  fact  sensitive  to  X  .  In  practice  a  great 
deal  of  experimental  information  is  available  (such  as  the  time 
distribution  at  each  value  of  c)  and  the  parameters  can  be 
decoupled  and  determined  independently.  An  additional  important 
piece  of  information  that  can  be  extracted  is  the  concentration 
of  deuterium  originally  present  in  the  hydrogen;  on  the  basis 
of  gamma  ray  yield  this  is  determined  to  be  less  than  one  part 
per  million. 

The  value  of  X  which  is  obtained  from  analysis  of  the 
PP 

results  of  the  fusion  experiment1 is  1.89  x  106  ±0.4  sec-1. 

This  corresponds  to  a  partial  lifetime  of  the  p,p  atom  of  0.53 
p.sec.  As  has  been  pointed  out  earlier,  the  capture  rate  of 
muons  in  the  singlet  |xp  system  is  only  about  14%  higher  than 
the  capture  rate  when  they  are  bound  in  the  ortho- (p|j,p)  molecu¬ 
lar  ion;  this,  together  with  the  fact  that  in  the  present 
experiment  useful  neutron  data  can  be  obtained  only  after 
about  1.2  p,sec  following  the  arrival  time  of  the  muon,  implies 
that  the  muon  capture  rate  as  measured  in  the  present  experi¬ 
ment  is  quite  insensitive  to  the  value  of  X  .  Experiments 

Xrh* 

using  bubble  chamber  techniques  record  events  occurring  at 
any  time  with  respect  to  the  muon  stopping  time,  and  are  there¬ 
fore  somewhat  more  dependent  on  a  knowledge  of  this  transfer 


rate . 
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IV.  DESCRIPTION  OF  THE  EXPERIMENT 
A.  The  Meson  Beam  and  the  Geometrical  Arrangement 

A  magnetic  field-free  channel  placed  in  the  cyclotron 
fringe  field  permitted  the  extraction  of  a  lower  energy  muon 
beam  than  would  otherwise  be  possible.  A  128  MeV/c  muon  beam 
was  bent  through  41°  by  a  wedge  shaped  dipole  magnet.  The 
beam  then  passed  through  a  pair  of  quadrupole  magnets  which 
were  set  so  as  to  maximize  the  counting  rate  in  the  beam 
counter  telescope;  this  condition  corresponds  to  bringing  the 
beam  to  a  focus  within  the  hydrogen  target  (see  Fig.  5) . 

Beryllium  moderator  was  used  to  reduce  the  muon  energy 
from  the  initial  60  MeV  to  the  appropriate  energy  which  would 
allow  muons  to  stop  between  the  neutron  detectors  at  the  center 
of  the  target  (see  Figs.  4,  8).  Beryllium  was  chosen  as  the 
optimum  material  for  reducing  multiple  scattering  of  the  beam. 

It  was  required  that  muons  of  26  MeV  impinge  on  the  window  of 
the  target ;  the  24 . 5  MeV  muons  then  entering  the  hydrogen  would 
come  to  rest  at  the  proper  position.  The  differential  stopping 
distribution  of  the  muons  was  measured  as  a  function  of  the 
amount  of  moderator  placed  in  the  beam.  The  approximately 
correct  amount  of  moderator  was  chosen  and  was  then  varied 
slightly  while  muon  decay  electrons  were  counted;  the  final 
working  conditions  were  chosen  to  be  those  which  maximized 
the  countring  rate  of  decay  electrons  in  the  electron  telescope. 
Differential  and  integral  stopping  distributions  as  well  as 
the  decay  electron  counting  rate  are  shown  as  a  function  of 
moderator  thickness  in  Figs.  9  and  10. 
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The  differential  range  distribution  of  the  muon  beam  as 
it  appears  in  Fig.  9  has  a  full  width  at  half  maximum  of 
3.15  grams  of  carbon.  This  width  is  due  to  the  momentum  selec¬ 
tion  properties  of  the  magnet  and  beam  channel;  interpreted  in 
this  way  the  beam  is  seen  to  have  a  momentum  dispersion  of 
±  5  MeV/c  or  ±  3.8%.  At  the  energies  of  interest  hydrogen 
has  a  stopping  power  2.35  times  greater  than  has  carbon.  The 
3 -gram  full  width  as  measured  with  carbon  then  translates  into 
7.3  in.  of  liquid  hydrogen.  The  neutron  detectors  are  each 
6  in.  long  in  the  direction  of  the  target  axis  so  that  most  of 
the  muons  are  expected  to  stop  directly  in  front  of  them. 

Beam  distribution  studies  in  a  styrofoam  dummy  target  which 
simulates  the  stopping  power  of  liquid  hydrogen  confirm  the 
expected  stopping  pattern  (see  Appendix  II) .  The  range  distri¬ 
bution  as  measured  with  decay  electrons  is ,  of  course ,  broader 
since  the  detectors  subtend  a  large  solid  angle  even  when  the 
beam  stopping  point  has  shifted  slightly;  as  a  result  the 
counting  rate  of  electrons  and  neutrons  is  less  sensitive  to 
small  changes  in  beam  momentum  than  might  otherwise  be  expected. 
On  the  other  hand  small  momentum  shifts  in  the  beam  are  likely 
to  be  accompanied  by  changes  in  beam  purity  and  hence  in  the 
accidental  events  rate.  For  this  reason  the  bending  magnet 
current  was  continually  monitored  on  a  chart  recorder  during 
the  run.  The  apparent  shift  between  the  peak  of  the  electron 
curve  and  the  half  height  point  of  the  integral  stopping  curve 
in  Fig.  10  represents  the  thickness  of  the  hydrogen  between  the 
beam  counters  #3  and  the  projected  center  of  the  decay  electron 
telescope . 
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Two  different  beam  collimator  arrangements  were  used 
during  two  successive  weeks  of  running  the  muon  capture  experi¬ 
ment.  One  was  characterized  by  an  aperture  4  in.  in  diameter, 
the  other  by  a  3  in.  opening;  although  the  latter  reduced  the 
number  of  muons  entering  the  target,  a  somewhat  improved  ratio 
of  wall  stopping  muons  to  hydrogen  stopping  muons  was  observed. 
The  data  from  these  two  independent  runs  were  analyzed  separately 
and  showed  no  significant  difference;  this  is  strong  confirmation 
that  correct  account  was  made  of  phenomena  relating  to  wall 
stoppings  and  wall  penetration.  The  two  geometrical  arrange¬ 
ments  are  shown  in  Figs.  8a  and  8b;  they  will  be  referred  to 
as  the  "3  inch"  and  the  "4  inch"  conditions,  respectively. 

The  meson  beam  telescope  normally  consisted  of  four  plastic 
scintillation  detectors  each  viewed  by  a  photomultiplier.  Two 
extra  scintillation  counters  were  inserted  for  special  purposes; 
one  to  serve  as  the  initial  counter  in  a  time -of -flight  system, 
and  the  other  to  be  in  anti-coincidence  with  the  beam  telescope 
for  differential  range  studies.  Table  III  lists  the  detectors 
along  with  some  of  their  properties . 

Figures  5  and  8  show  the  locations  of  the  beam  detectors 
with  respect  to  the  shielding  arrangement.  Counter  #1  was 
embedded  in  the  front  shielding  wall  but  could  be  withdrawn 
from  above  to  check  line-up.  It  was  thin  to  prevent  beam 
scattering,  a  problem  which  rendered  it  impractical  to  employ 
the  time-of-flight  system  as  an  added  requirement  during  actual 
running.  When  the  purified  muon  beam  is  used  the  pions  are  of 
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lower  momentum  than  the  muons  which  are  accepted  by  the  colli¬ 
mator  system.  They  are,  therefore,  bent  through  a  larger  angle 
and  most  of  them  strike  the  shielding  blocks  or  the  sides  of 
collimator  walls.  The  location  of  counter  #1  was  chosen  so  as 
to  intercept  those  pions  which  were  likely  to  give  rise  to  a 
detected  neutron.  These  events  were  discarded,  as  shall  be 
described,  reducing  the  accidental  neutron  rate  quite  substan¬ 
tially. 

Counters  #3a  and  #3b  were  located  adjacent  to  one  another 
downstream  of  the  moderator;  the  second  of  the  two  was  extremely 
thin  (0.02  in.)  to  reduce  the  number  of  muons  stopping  in 
plastic  to  a  tolerable  level.  The  two  counters  were  placed 
in  coincidence  with  one  another  so  as  to  reduce  the  accidental 
counting  rate  in  the  beam  telescope.  An  accidental  coincidence 
in  this  case  would  very  likely  be  due  to  a  muon  stopping  in  the 
beryllium  moderator  coincident  with  a  noise  pulse  and  would  be 
electronically  indistinguishable  from  a  hydrogen  stopping. 

Counter  #3b  consisted  of  a  20  mil  sheet  of  plastic  scintillator 
in  a  protective  frame  which  supported  a  loose  envelope  of  thin 
aluminum  foil  (to  provide  for  light  collection)  and  a  1  mil 
sheet  of  black  mylar. 

B.  Electron  and  Neutron  Detectors 

The  hydrogen  target  was  surrounded  by  four  11|"  x  14"  x  3/8" 
plastic  scintillation  counters;  behind  each  of  these  was  mounted 
a  neutron  detector  surrounded  by  a  least  2  in.  of  lead  on  five 
of  its  sides.  The  thin  counters  were  used  in  coincidence  with 
the  neutron  detectors  to  define  electrons  and  in  anti-coincidence 
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to  select  neutral  particles.  These  counters  were  labeled  A^, 

A2,  A3,  and  A^  in  a  clockwise  direction  beginning  at  the  upper 
right  facing  the  beam;  the  respective  neutron  counters  were 
called  .  The  anti-coincidence  counters  were  constructed  so 
as  to  form  four  sides  of  a  box  around  the  target;  they  subtended 
approximately  70%  of  the  solid  angle  available  to  decay  electrons 
from  stopping  muons.  On  the  side  facing  the  target  each  counter 
was  covered  only  with  0.00075  in.  of  aluminum  foil  and  0.001  in. 
black  mylar;  these  precautions  were  taken  lest  muons  stop  in 
low  atomic  number  materials  and  be  captured  there.  The  anti- 
coincidence  counters  were  used  in  four  ways  in  the  experimental 
logic:  1.  As  a  prompt  anti-coincidence  at  muon  arrival  time 

to  warn  of  a  muon  having  left  the  target  and  having  penetrated 
the  scintillator;  2.  In  coincidence  with  the  neutron  detectors 
to  define  electrons;  3.  In  prompt  anti-coincidence  with  a  pulse 
in  the  neutron  detector  to  exclude  electrons  from  consideration 
as  events;  4.  As  a  10  p,sec  anti-coincidence  (on  film)  to  reduce 
accidental  background. 

The  neutron  detectors  were  made  of  liquid  scintillator 
NE213  which  was  purified,  deoxygenated  and  sealed  in  specially 
constructed  glass  vessels  6"  x  5"  x  4"  in  size.  An  EMI  photo¬ 
multiplier,  9530B,  directly  viewed  a  5"  x  6"  face  of  the 
rectangular  volume.  The  1.8  liter  container  was  immersed  in 
magnesium  oxide  powder  inside  a  thin  iron  box;  each  of  the  four 
assembled  counters  was  inserted  into  a  lead  well  around  the 
hydrogen  target.  The  largest  single  source  of  counts  in  the 
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neutron  counters  are  gamma  rays  due  to  bremsstrahlung  of  the 

muon  decay  electrons  in  their  passage  through  the  high  Z  target 

walls.  In  order  to  distinguish  between  neutron  and  gamma  ray 

induced  recoils  in  the  detector,  pulse  shape  analysis  of  the 

25 

form  developed  by  Brooks  was  employed.  The  light  produced 

by  a  charged  particle  in  a  scintillator  is  composed  of  a  long 

lived  (~  300  nsec)  component  superposed  on  the  more  familiar 

10  nsec  initial  pulse.  The  size  of  the  slower  component  is 

always  proportional  to  energy  loss,  while  the  fast  pulse  is 

linear  only  for  minimum  ionizing  particles  and  saturates  for 

dE 

particles  with  a  high  rate  of  energy  loss  ^  such  as  1-5  MeV 
protons.  Thus,  one  can  recognize  neutrons  and  gamma  rays  by 
the  ratio  of  the  two  components  of  the  light  pulse  produced  by 
the  proton  recoil  or  Compton  recoil.  Appendix  I  describes  the 
saturation  effect  which  enables  particle  discrimination  and 
which  also  causes  the  non-linear  energy  response  in  the  case 
of  protons.  The  circuit  in  Fig.  11  was  used  to  provide  a 
particle  discrimination  pulse,  derived  from  the  slow  component, 
which  gated  neutral  events  and  was  also  displayed  on  an  oscil¬ 
loscope  to  be  photographed  and  interpreted  visually.  Two 
pulses  of  opposite  polarity  are  taken  from  the  anode  and  dynode 
of  the  photomultiplier;  the  positive  pulse  enters  a  low  pass 
filter  to  extract  the  fast  component,  it  then  discharges  a 
capacitor  which  was  charged  to  a  negative  voltage  by  the  fast 
pulse  from  the  anode.  The  system  is  balanced  so  that  in  the 
case  of  gamma  rays  there  is  very  little  residual  voltage  on 
the  capacitor,  while  in  the  case  of  neutrons  there  is  net 
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positive  charge .  The  voltage  at  the  capacitor  is  sensed  by  an 
emitter-follower  and  is  then  amplified  and  sent  to  a  discriminator 
This  residue  pulse  can  be  displayed  as  a  function  of  particle 
energy  by  applying  it  and  the  fast  component  pulse  to  the 
horizontal  and  vertical  deflection  plates  of  a  modified  oscil¬ 
loscope  whose  electron  beam  is  unblanked  only  when  the  chosen 
logical  requirements  have  been  met  (see  Fig.  12) .  This  sytem 
of  display  is  useful  in  adjusting  the  pulse  shape  balance  con¬ 
trols  and  the  residue  discriminator  and  for  periodic  checking 
of  the  particle  discrimination  quality. 

The  pulse  shape  discrimination  circuit  shown  in  Fig.  11 
was  located  directly  at  the  base  of  each  neutron  detector  photo¬ 
multiplier.  The  output  amplification  stage  was  not  only  neces¬ 
sary  to  provide  enough  gain  to  send  the  signal  up  to  the  count¬ 
ing  area,  but  was  biased  so  as  to  saturate  on  the  unwanted 
initial  negative  spike  and  suppress  it  relative  to  the  somewhat 
delayed  neutron  residue.  The  photomultiplier  anode  and  dynode  11 
supplied  pulses  to  this  circuit;  dynodes  9  and  10  supplied  fast 
pulses  (actually  the  tube  has  a  25  nsec  rise  time)  one  of 
which  went  directly  to  the  oscilloscopes  for  photographing  and 
the  other  was  further  amplified  and  mixed  with  its  sister  pulses 
from  the  other  three  detectors.  This  mixed  pulse  then  went  to 
two  pulse  height  discriminators;  one  served  as  a  lower  threshold, 
the  other  as  an  upper  cutoff  on  neutron  events  (but  not  elec¬ 
trons)  .  The  upper  discriminator  which  was  set  for  about  15  MeV 
neutrons,  eliminated  high  energy  bremsstrahlung  events  whose 
slow  residue  might  be  sufficiently  large  to  pass  the  residue 
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discriminator  and  be  counted  as  neutrons.  Discriminator 
settings  and  photomultiplier  voltage  stability  were  periodi¬ 
cally  checked  by  displaying  the  oscilloscope  input  pulse  on  a 
kicksorter  which  was  gated  by  the  logic  output  (discriminators, 
etc.) .  Such  pulse  height  spectra  (Fig.  24a)  were  also  used 
to  measure  the  intrinsic  energy  resolution  of  the  detectors. 
With  low  energy  gamma  sources  simulating  proton  recoils  in  a 
detector  of  the  volume  used  (4"  x  5"  x  6")  the  resolution 
spread  is  almost  entirely  due  to  statistical  fluctuations  in 
the  number  of  photoelectrons  emitted  at  the  tube  photocathode. 
The  resolution  was  taken  to  be  ±  10%  at  2.76  MeV  (gamma  ray 
energy)  and  proportional  to  1//E.  Figure  12  shows  the  energy 
versus  residue  plot  of  a  Pu-Be  source,  showing  separated 
neutron  and  gamma  ray  bands.  Gamma  ray  sources  are  used  to  set 
the  energy  scale  which  is  nonlinear  for  proton  recoils.  Cali¬ 
bration  of  the  scale  involved  an  independent  experiment  which 
is  described  in  Appendix  I.  Neutron-gamma  ray  separation  is 
possible  only  above  some  minimum  energy  depending  on  the  light 
collection  efficiency  of  the  detector  and  other  factors.  Fluc¬ 
tuations  in  the  number  of  photoelectrons  received  from  the  slow 
component,  which  contain  less  than  one-third  the  total  output 
charge,  result  in  an  inability  to  resolve  the  two  bands  when 
the  light  level  is  low. 

C.  Logic  and  Data  Storage 

Three  types  of  coincidence  circuits  were  used  in  the 
electronic  logic  of  the  muon  capture  experiment.  Fast  coinci¬ 
dences  between  pairs  of  beam  counters  were  made  in  circuits  of 
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the  Franzini  type  adjusted  to  have  resolving  times  of  about 

10  nsec.  The  use  of  the  improved  duty-cycle  and  purified  beam 

eliminated  the  need  for  ultra-fast  circuitry  while  the  duration 

of  the  run  made  it  convenient  to  use  broader  resolving  times 

which  could  tolerate  small  drifts  in  photomultiplier  voltages, 

27 

etc.  Circuits  of  the  type  described  by  Fitch  were  used  with 

30-60  nsec  resolving  times  for  making  anticoincidences  and 

coincidences  between  outputs  of  other  circuits .  The  slow  p,sec 

logic  needed  for  making  delayed  coincidences  between  muon  pulses 

and  electrons  and  neutrons  employed  a  direct  coupled  system  of 

modular  gates  and  univibrators  developed  for  general  use  in  the 
28 

laboratory.  Photomultiplier  pulses  were  led  from  the  detectors 
in  the  cyclotron  experimental  area  to  the  laboratory  counting 
area  on  125  ohm  (RG-63U)  cable?  pulses  to  the  time -of -flight 
system,  however,  were  sent  via  low  loss  styrofoam  filled  cable. 
Intercircuit  connections  were  usually  made  with  50  ohm  cable 
(RG-58)  while  long  delays  were  achieved  by  reshaping  the  output 
of  1500  ohm  delay  line.  Fast  pulses  destined  for  oscilloscope 
photography  were  delayed  with  RG-63U.  Transistorized  circuitry 
was  used  throughout  except  for  high  voltage  power  supplies  and 
Hewlett-Packard  distributed  amplifiers  used  in  mixing  fast 
pulses.  Figure  14  is  a  simplified  block  diagram  of  the  logic. 

The  following  eleven  logical  combinations  of  detector  outputs 
were  formed  and  recorded  on  a  set  of  scalers: 
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12 

123a3b 

123a3b(EAi) 

(ICN^  (EA±) 

(SN.)e 
1  u 

(non-e) (Res) 

(mu)  e 

u 

(mu)  n 


a  gated  train  of  oscillator  pulses 
a  coincidence  between  first  two  beam 
counters 
the  muon  beam 

muons  not  leaving  target  ="mu  stop"  =  |i 

decay  electrons  =  e^ 

neutral  particles  converting  in  the 

neutron  detector  =  non-e 

neutrons,  (slow  component  residue  is 

sufficiently  large)  =  n 

muon  gated  electrons,  real  and  accidental 

=  e  ,  e  (e  -  e  =  e) 
r  a  r  a 

gated  neutrons  =  n  ,  n  =  "events " . 
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The  notation  (£A^)  or  (EN^)  refers  to  the  mixture  of  the  four 
detector  outputs.  In  addition  to  the  requirements  noted  above, 
each  "|i"  signal  had  to  pass  the  "beam  gate";  that  is,  it  was 
used  only  if  it  occurred  during  the  8  millisec  when  the  meson 
beam  was  being  ejected  with  good  duty-cycle.  The  requirement, 
nl,  discarding  events  which  are  in  coincidence  with  an  incoming 
pion,  was  omitted  from  the  logic  but  was  enforced  when  the 
scanned  oscilloscope  photographs  were  analyzed.  As  will  be 
seen  later,  there  was  often  a  delay  between  the  detection  of  an 
incoming  pion  and  the  detection  of  the  neutron  due  to  the  flight 
time  of  the  neutron  from  the  point  at  which  it  was  produced 
upstream  to  the  detector  within  the  shielding  house.  A  5  MeV 
neutron  takes  32  nsec  to  traverse  a  meter;  so  that  while  the 
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flight  time  is  negligible  for  neutrons  originating  in  the  hydro¬ 
gen  target,  it  is  not  when  the  neutron  is  produced  in  the  walls 
of  a  beam  collimator  early  in  the  beam  path. 

To  permit  careful  time  and  energy  measurements  as  well  as 
to  check  neutron-gamma  ray  discrimination  and  anticoincidence 
operation,  several  detector  pulses  were  photographed.  A  sequence 
of  mirrors  was  used  to  record  simultaneously  the  traces  of  the 
Tektronix  oscilloscopes  517  and  555  on  film.  Both  the  oscil¬ 
loscopes  and  the  film  advance  mechanism  were  triggered  on  the 
arrival  of  either  an  n  or  n  signal;  that  is,  a  muon  gated 
neutron  that  occurred  during  a  5  [isec  interval  beginning  either 
1/2  |j,sec  after  a  muon  pulse  or  15  p.sec  thereafter.  The  former 
interval  includes  events  from  muon  capture  as  well  as  some 
random  background;  the  latter  interval  only  yields  background 
events  since  all  captures  and  decays  are  assumed  to  have  taken 
place  before  15  p,sec  have  elapsed.  The  background  as  measured 
in  the  second  interval  can  be  used  to  correct  the  number  of 
events  measured  in  the  first  interval,  as  the  background  is 
uncorrelated  in  time  with  individual  muons. 

Figure  15  shows  typical  oscilloscope  traces.  The  "517" 
oscilloscope  trace  records  the  pulse  from  the  neutron  detector 
and  a  pulse  from  beam  counter  #1;  these  are  used  to  determine 
whether  a  pion  in  the  beam  might  have  given  rise  to  the  neutron 
event.  The  other  oscilloscope  displays  two  traces,  of  which 
the  lower  is  a  record  of  the  neutron  detector  output  (the  fast 
pulse  amplitude  and  the  slow  residue)  used  to  measure  the 
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neutron  energy  and  to  perform  neutron-gamma  ray  discrimination. 
The  upper  trace  carries  the  neutron  pulse  again  in  attenuated 
form  and  a  shaped  "p"  pulse;  it  may  also  display  pulses  in  the 
anticoincidence  counters  if  there  should  be  any.  Near  the  end 
of  the  trace  is  a  pulse  whose  length  identifies  one  or  another 
of  the  four  neutron  detectors  as  the  source  of  the  event. 

Simultaneously  with  data  being  taken  on  film,  a  "digitron" 
system  received  pulses  representing  the  muon  and  electron  and 
stored  their  time  interval  distribution  in  the  magnetic  core 
memory  of  a  400  channel  pulse  height  analyzer,  making  possible 
lifetime  measurements  of  very  high  statistical  accuracy. 

D.  Data  Taking  Procedure 

Because  of  the  low  muon  capture  rate  as  compared  with 
decay,  only  10-20  events  were  expected  per  hour,  necessitating 
a  run  of  two  weeks  to  collect  sufficient  data  and  to  obtain 
accurate  knowledge  of  the  corrections.  The  "4  inch"  geometrical 
arrangement  (Fig.  8b)  was  used  during  the  first  week  of  running; 
under  these  conditions  400  muons  per  second  stopped  in  hydrogen. 
Under  "3  inch"  (Fig.  8a)  conditions  used  during  the  second  week 
only  280  muons  per  second  could  be  stopped  in  hydrogen.  Fre¬ 
quent  checks  were  made  to  ensure  the  constancy  of  certain  ratios 
of  counting  rates,  such  as  e/p.,  12/Time,  p/12,  ea/er»  etc.  These 
ratios  measure  the  beam  purity,  changes  in  the  muon  stopping 
distribution,  vibrating  target  efficiency,  and  general  stability 
of  the  apparatus.  The  response  of  the  neutron  detectors  was 

periodically  checked  with  a  series  of  gamma  ray  sources  whose 

.  137 

energy  was  m  the  region  of  interest;  Cs 


-0.67  MeV, 
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Co^  -  1.25  MeV  (the  detectors  are  not  able  to  resolve  the 

24 

two  closely  spaced  lines),  Na  -  1.35  MeV  and  2.76  MeV.  The 
neutron-gamma  ray  discrimination  quality  was  checked  every  few 
hours  using  a  plutonium-beryllium  source  which  emits  neutrons 
of  energy  up  to  about  6  MeV  and  gamma  rays  of  4.4  MeV. 

It  was  realized  during  an  earlier  version  of  the  muon  cap¬ 
ture  experiment  that  a  significant  source  of  background  was 
produced  by  decay  electron  bremsstrahlung.  Since  the  effect 
was  electron  associated,  it  was  possible  to  measure  it  quite 
accurately  by  stopping  positive  muons  in  the  target  and  look¬ 
ing  for  neutrons  under  exactly  the  same  conditions  as  with  the 
negative  beam.  It  is  rather  easy  to  reverse  the  beam  polarity, 
one  has  only  to  reverse  the  cyclotron  magnetic  field  direction 
so  that  mesons  of  the  opposite  sign  are  extracted,  the  bending 
magnet  and  quadrupole  magnet  polarities  are  also  reversed,  but 
nothing  in  the  geometrical  arrangement  has  to  be  changed.  Fre¬ 
quent  alternations  of  beam  polarity  also  made  possible  a  very 
precise  measurement  of  the  positive  and  negative  muon  lifetime 
ratio  confirming  the  prediction  of  particle-antiparticle  life¬ 
time  equality  by  the  TCP  theorem.  This  will  be  the  subject  of 
a  separate  report. 

Whenever  data  were  being  accumulated  the  count  scalers  and 
the  oscilloscope  system  were  stopped  once  pa: hour?  the  numbers 
of  counts  that  were  registered  in  each  of  the  eleven  channels 
were  recorded  and  the  scalers  were  reset  to  zero;  a  record  was 
made  of  the  film  frame  number  so  that  a  strict  correlation 
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could  be  established  between  the  photographed  events  and  the 
scaler  record.  Three  times  each  day  radioactive  sources  were 
installed  and  the  oscilloscope  "events"  trigger  was  simulated 
so  that  photographs  could  be  taken  of  the  neutron  counter  response 
to  the  various  sources.  Pictures  were  also  taken  of  decay  elec¬ 
tron  events  by  relaxing  the  anticoincidence  requirement;  these 
were  useful  in  establishing  the  pulse  size  due  to  minimum 
ionizing  particles  in  the  anticoincidence  counters. 

V.  DATA  ANALYSIS 

A.  Film  Scanning  and  Data  Handling 
During  the  course  of  the  muon  capture  run,  about  6000  ft 
of  35  mm  film  were  used  in  photography  of  the  oscilloscope 
traces.  This  included  data  taken  with  the  positive  and  nega¬ 
tive  beam  as  well  as  radioactive  source  calibration  pictures. 

The  film  was  reprojected;  pulse  sizes  and  spacings  were  measured 
using  a  template  and  the  measurements  were  punched  on  IBM  cards, 
one  card  per  frame.  The  following  information  was  recorded: 
the  film  frame  number,  a  code  denoting  the  type  of  data,  the 
number  of  the  neutron  detector  responsible  for  the  event,  the 
size  of  the  fast  output  pulse  from  the  neutron  detector,  the 
size  of  the  fast  pulse  which  was  attenuated  by  a  factor  of  4 
(to  increase  the  dynamic  range  of  the  oscilloscope),  the  size 
of  the  slow  component  residue,  the  time  between  the  neutron 
pulse  and  the  muon  trigger,  the  time  between  the  neutron  pulse 
and  anticoincidence  counter  pulse  (if  any),  the  time  between 
the  neutron  pulse  and  the  pulse  from  beam  counter  #1  (if  any)  , 
and  a  code  denoting  whether  the  event  occurred  during  the  "real " 
gating  interval  or  during  the  delayed  "accidental"  interval. 


-39- 


All  of  the  above  information  except  the  counter  #1  pulse  and 
its  time  marker  appeared  on  the  two  traces  of  the  Tektronix  555 
oscilloscope  with  a  sweep  speed  of  1  jj,sec  per  centimeter  (10  p,sec 
full  scale)  ;  the  517  oscilloscope  displayed  1  |j.sec  of  time  in 
the  vicinity  of  the  neutron  pulse .  From  time  to  time  segments 
of  the  film  were  rescanned  to  check  scanning  efficiency  and 
accuracy;  absolute  time  calibration  was  made  possible  by  photo¬ 
graphing  a  crystal  controlled  oscillator  output  and  pulse  trains 
of  various  amplitudes . 

The  data  cards  were  sorted  according  to  neutron  detector 
number  and  the  presence  or  absence  of  a  "pion"  pulse  (counter  #1) 
at  the  time  of  the  neutron  (or  shortly  before) .  The  proper  time 
interval  to  exclude  was  chosen  on  the  basis  of  a  distribution 
curve  of  pi on-neutron  time  intervals  plotted  from  the  data. 

This  plot  also  indicates  the  fraction  of  events  that  were  need¬ 
lessly  excluded  because  of  a  random  pulse  in  counter  #1;  these 
accidentals  are  taken  into  consideration  in  the  analysis.  The 
fast  trace  also  pictures  real  events  which  occurred  near  the 
beginning  of  the  gating  interval  and  can  be  used  to  obtain  a 
high  resolution  time  interval  distribution  of  the  early  neutrons 
which  are  largely  from  captures  in  the  iron  walls  of  the  target. 
Since  the  gating  interval  does  not  begin  until  about  1/2  pisec 
after  the  muon  stopping  time,  the  captures  in  the  very  high  Z 
materials  such  as  zirconium  and  silver  are  not  visible  in  the 
oscilloscope  data;  these  capture  events  are  clearly  observed, 
however,  in  neutron  time  distribution  data  taken  with  the  digi- 
tron.  Early  time  data  is  used  to  obtain  an  accurate  estimate 
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of  the  number  of  muons  stopping  in  the  target  walls;  this  number 
while  not  interesting  in  itself  is  important  in  deducing  the 
number  of  muons  which  penetrate  the  walls  and  stop  in  the 
lower  Z  coating  of  the  anticoincidence  detectors . 

Events  which  contain  evidence  of  a  pulse  in  the  anticoinci¬ 
dence  counter  are  excluded  for  several  reasons?  a  fraction  of 
these  events  appeared  only  because  the  NA  anticoincidence 
failed  due  to  the  small  size  of  the  pulse  in  the  A  counter,  a 
fraction  appeared  because  the  (mu) A  anticoincidence  failed  for 
the  same  reason,  and  finally,  some  events  appeared  because  of 
an  accidental  (random)  neutron  pulse  in  conjunction  with  the 
actual  decay  electron  of  an  actual  muon.  Thus  the  visual  exclu¬ 
sion  of  events  with  A-pulses  not  only  increases  the  efficiency 
of  the  electronic  anticoincidences,  but  also  eliminates  a  sub¬ 
stantial  fraction  of  random  events  with  the  aid  of  the  very 
large  solid  angle  of  the  A  counters  and  the  10  p.sec  of  time 
which  are  visible  on  the  film.  This  requirement  is  analogous 
to  the  condition  in  a  bubble  chamber  experiment  that,  if  a  muon 
is  to  be  captured,  there  should  at  least  be  no  decay  electron 
track  emerging  from  its  stopping  point. 

For  each  individual  neutron  detector,  for  the  data  accumu¬ 
lated  during  each  day  of  running,  a  plot  was  made  of  the  size 
of  the  fast  pulse  in  the  neutron  detector  against  the  size  of 
the  slow  residue.  This  plot  is  analogous  to  the  type  of  oscil¬ 
loscope  display  that  was  available  during  the  run.  An  example 
of  such  a  plot  along  with  similar  plots  for  radioactive  sources 
is  given  in  Fig.  12c,  d.  By  studying  these  plots  one  could 
eliminate  all  events  which  fell  in  the  region  of  gamma  rays? 
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one  could  also  determine,  for  each  detector,  the  lowest  energy 
that  would  permit  reliable  neutron-gamma  ray  discrimination. 

Events  which  passed  the  tests  mentioned  thus  far  could  no  longer 
be  rejected  on  an  individual  basis;  each  one  had  the  charac¬ 
teristics  of  muon  capture  in  hydrogen.  These  events  were  stored 
on  magnetic  tape  for  subsequent  handling  on  an  IBM  1401  Data 
Processing  System,  they  were  still  grouped  according  to  neutron 
detector  number  and  day  of  running;  data  from  the  3  and  4  inch 
geometrical  situations  were  kept  separate  throughout.  Muon- 
neutron  time  distributions  were  plotted  from  these  events.  These 
are  shown  for  the  two  geometries  in  Pigs.  16a,  16b.  They  are 
characterized  by  an  approximately  2.2  |j.sec  lifetime  component, 
as  is  expected  from  muon  capture  in  hydrogen,  and  a  0.2  p.sec 
component  from  iron  captures.  Accidental  events  whose  time 
distribution  is  shown  in  Fig.  17  were  subtracted  from  the  life¬ 
time  plot.  For  both  the  3  inch  and  4  inch  running  conditions 
a  time  delay  was  chosen  (using  the  lifetime  plots)  beyond  which 
the  effect  of  iron  captures  would  be  small  but  not  necessarily 

negligible  (since  a  sufficiently  accurate  correction  is  possible) , 

.  .  .  2 

The  time  distributions  were  fitted,  by  minimizing  y  ,  by  a 
function  of  essentially  the  form  exp(-t/0.2)  +  Ng  exp(-t/x) 
with  N^,  N2  and  x  as  parameters.  Actually  the  function  was 
modified  slightly  to  take  into  account  the  finite  resolution  of 
the  grouped  data.  The  results  of  this  fit  are: 

x(Lisec) 

3  Inch  2.13  ±  0.15 


4  Inch 


2.10  ±  0.12 
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The  data  from  positive  beam  running  were  treated  similarly, 
the  time  distributions  are  shown  in  Fig.  18.  There  is,  of 
course,  no  short-lived  component  in  the  lifetime  plot;  the 
results  are: 

x(Lisec) 

3  Inch  2.24  ±  0.17 

4  Inch  2.21  ±  0.15 


B.  Calculation  of  Capture  Rate 
The  experimental  value  of  the  muon  capture  rate  is  given 


by 


Rate  = 


N-c 


M  ‘  Xo 


E*SA*T 

where  N  is  the  number  of  detected  neutrons ,  c  represents  neu¬ 
trons  due  to  sources  other  than  hydrogen,  E  is  detector  effi¬ 
ciency,  SA  is  detector  solid  angle,  and  M  is  the  number  of 
muons  stopping  in  hydrogen,  T  is  the  time  gating  factor.  Each 
of  these  factors  will  be  considered  in  turn. 

1.  Muons  Stopping  in  Hydrogen 
One  can  arrive  at  the  number  of  muons  stopping  in  hydrogen 
from  two  directions;  the  number  of  (123a3b)A  counts  or  the 
number  of  (e  -  e  )  decay  electrons.  Several  corrections  must 
be  made  to  the  raw  number  of  beam  counts  to  obtain  a  muon  stop¬ 
ping  number;  these  are  due  to  muons  stopping  in  the  front  wall 
of  the  target,  the  side  walls  and  in  the  collimator  before  the 
target.  Experimental  runs  with  a  dummy  target  (see  Appendix  II) 
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were  carried  out  in  order  to  measure  the  muon  stopping  distribu¬ 
tion  and  to  determine  the  percentage  of  "p,"  counts  which  are 
muons  entering  the  target  region.  It  was  found  that  only  about 
three-quarters  of  the  (123a3b)A  triggers  represent  muons  which 
leave  the  last  collimator,  the  rest  having  trajectories  which 
failed  to  enter  the  target.  Because  of  this  large  loss  factor, 
this  method  can  only  serve  as  a  rough  check  and  it  was  decided 
to  determine  the  number  of  muons  stopping  in  hydrogen  solely 
from  the  number  of  detected  decay  electrons .  This  second  method 
is  more  reliable  in  any  case .  The  important  advantages  are : 
i)  except  for  a  small  and  known  number  of  accidentals  a  (mu)e 
coincidence  can  only  be  a  decay  electron;  ii)  electrons  are 
required  to  pasb  the  same  time  gate  as  neutrons;  iii)  electrons 
are  detected  with  the  same  neutron  detectors  (and  hence  almost 
the  same  solid  angle)  as  are  neutrons;  iv)  since  the  time  gate 
is  sufficiently  delayed  only  decay  electrons  from  muons  stopping 
in  hydrogen  are  counted.  The  number  of  muons  stopping  in  hydro¬ 
gen  is  given  by 

M  =  (e  -  eJ/T-SA-L 

i  a 

where  T  is  the  time  gate  factor,  SA  is  the  solid  angle  for  elec¬ 
tron  detection,  and  L  is  a  loss  factor  for  electrons  due  to 
target  side  walls . 

Table  IV  lists  the  various  factors;  (e  -  e  )  has  been  set 

r  a 

equal  to  e .  In  the  negative  beam  e  /e  is  1.5%  and  2.1%, respec- 
tively,  for  the  3  and  4  inch  geometries.  The  time  factor  was 
deduced  from  the  lifetime  plot  (Fig.  16)  assuming  a  decay  life¬ 
time  of  2.2  p,sec  and  remembering  that  the  same  muon-generated 
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univibrator  output  gates  both  electrons  and  neutrons  (Fig.  14) . 
The  solid  angle  for  detection  of  an  electron  assuming  the 
known  muon  stopping  distribution  in  hydrogen  (see  Appendix  II) 
was  calculated,  using  the  same  Monte  Carlo  code  as  was  used 
for  the  neutron  detector  efficiency  calculation  (see  Appendix  III) 
a  small  correction  was  applied  to  account  for  edge  effects  in 
the  neutron  detector. 

The  electron  loss  factor,  L,  consists  of  several  effects 
of  which  the  biggest  and  most  difficult  to  calculate  is  radia¬ 
tion  in  the  high  Z  target  walls.  A  fraction  of  the  decay  elec¬ 
trons  whose  spectrum  extends  up  to  53  MeV  will  radiate  suffi¬ 
cient  energy  that  the  remaining  electron  will  not  be  energetic 
enough  to  penetrate  the  neutron  detector,  and  will  therefore 
not  register  as  an  NA  coincidence.  This  effect  was  directly 
measured  with  the  dummy  target  and  simulated  target  walls; 
the  ratio  e/mu  was  recorded  with  the  walls  in  place  and  with 
the  walls  removed,  while  the  muons  maintained  the  same  stopping 
distribution  in  the  styrofoam  pseudohydrogen.  It  was  found 
that  5 %  of  the  electrons  were  lost  by  this  mechanism.  If  the 
energy  threshold  for  electrons  in  the  neutron  detector  is 
0.6  MeV  and  account  is  taken  of  the  intervening  materials  out¬ 
side  the  target,  then  a  plot  of  the  decay  electron  energy 
spectrum  indicates  that  another  5%  of  the  electrons  which  are 
headed  toward  the  detector  fail  to  be  counted.  Thus  a  loss 
factor  of  10%  ±  5%  is  used  to  correct  the  observed  number  of 


electrons. 
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2.  The  Time  Gate  Factor 

From  the  neutron  time  distribution  plot  a  time  interval 
was  chosen  which  excluded  most  of  the  iron  capture  events,  only 
the  data  within  this  new  time  interval  was  analyzed  to  extract 
the  muon  capture  rate .  Events  occurring  at  earlier  times  were 
isolated  and  used  to  provide  an  energy  spectrum  of  neutrons 
from  muon  capture  in  iron.  The  time  intervals  used  are  listed 
below: 


Start  of  Inter- 

End  of  Inter- 

Gate 

Geometry 

val  ( nsec) 

val  (usee) 

Factor  (T) 

3  Inch 

1.2 

5.2 

0.486  ±  1% 

4  Inch 

1.3 

5.1 

0.456  ±  1% 

3.  Solid  Angle  and  Detector  Efficiency 
The  total  solid  angle  for  the  four  neutron  detectors  was 
calculated  (Appendix  III)  to  be  19.6%  given  the  measured  muon 
stopping  distribution.  The  solid  angle  in  this  sense  refers  to 
the  probability  that  isotropically  directed  particles  emitted 
from  the  muon  stopping  point  will  enter  the  front  face  of  the 
neutron  detector.  In  the  case  of  neutrons  there  are  several 
other  effects  to  consider:  First,  the  neutron  may  be  transmitted 
through  the  lead  surrounding  the  detectors  and  enter  from  an 
unusual  direction,*  this  is  particularly  unlikely  because  the 
muon  stopping  distribution  is  peaked  in  front  of  the  detectors. 
Second,  a  neutron  has  a  substantial  probability  of  traversing 
the  entire  neutron  detector  without  inducing  a  proton  recoil; 
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the  detector  thickness  is  somewhat  less  than  one  mean  free  path 
for  a  collision  with  the  hydrogen  in  the  scintillator  volume  so 
that  some  neutrons  will  go  undetected.  Third,  even  if  a  colli¬ 
sion  with  a  proton  occurs,  the  recoil  proton  energy  may  be  too 
small  to  produce  a  pulse  that  can  pass  the  lower  energy  discrimi¬ 
nator  threshold. 

A  complete  calculation  of  neutron  detector  efficiency  was 
essential;  the  calculation  would  not  have  been  difficult  if  it 
were  not  for  several  complicating  conditions: 

i)  The  neutron  detectors  were  surrounded  by  lead  and  other 
materials  which  might  "reflect"  neutrons. 

ii)  The  path  length  in  the  detector  is  direction  dependent. 

iii)  Neutron  cross  sections  are  energy  dependent. 

iv)  The  detectors  are  thick  enough  to  permit  multiple  recoils. 

v)  The  hydrogen  in  the  target  itself  is  an  excellent  scatterer 
and  degrader  of  neutrons . 

vi)  The  response  of  the  detectors  is  non-linear  for  protons  of 
the  relevant  energies  (<5.2  MeV) . 

vii) The  detector  contains  carbon  as  well  as  hydrogen. 

For  these  reasons  a  complete  Monte  Carlo  calculation,  including 
all  of  the  above  effects,  was  necessary.  Appendix  III  describes 
the  calculation  which  yielded  a  distribution  of  the  number  of 
events  as  a  function  of  pulse  height.  In  the  limit  of  a  thin, 
isolated,  hydrogenous,  linear  detector  this  would  simply  be  the 
proton  single  recoil  energy  spectrum  which  would  in  fact  be  a 
flat  distribution  extending  from  zero  to  the  incident  neutron 


energy. 
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Since  the  calculation  gives  the  final  expected  yield  as  a 
single,  but  energy  cut-off  dependent,  quantity  it  is  not  mean¬ 
ingful  to  decouple  the  several  effects;  but,  because  the  elec¬ 
tron  yield  depends  on  the  solid  angle  also  (which  therefore 
"cancels  out"  in  the  final  capture  rate) ,  it  is  interesting  to 
list  the  extent  to  which  each  effect  determines  the  final 
neutron  yield.  The  numbers  are  approximate  and  are  not  needed 
in  the  analysis. 

Fraction  of  Neutrons  Observed  due  to 
_ Effect  of; _ 


Solid  angle 

0.20 

Hydrogen  degradation 

0.90 

Detector  mean  free  path 

0.50 

Energy  cut-off  ( 1  MeV) 

0.36 

The  corrected  observed  neutron  energy  spectrum  was  computed  as 
will  be  seen  below  and  the  calculated  distribution  was  fitted 
to  it,  normalizing  it  to  the  total  yield  above  a  cut-off 
energy  Eq.  The  efficiency  for  this  value  of  Eq  was  then  used 
to  obtain  the  capture  rate.  Figure  19  shows  the  corrected 
neutron  pulse  height  distributions  together  with  the  calculated 
pulse  height  spectrum  for  each  of  the  detectors  and  each  geome¬ 
trical  arrangement.  After  the  expected  spectrum  was  normalized 
2 

X  was  calculated  for  each  case.  The  errors  used  to  calculate 
2 

X  are  those  shown  on  the  graphs,  they  are  a  combination  of 
statistical  errors  and  propogation  errors  due  to  the  subtractions 
4.  The  Corrected  Number  of  Detected  Neutrons 
The  data  on  gamma  ray  sources  that  were  taken  during  each 
day  of  running  were  plotted  and  a  scale  factor  (to  convert  the 
arbitrary  units  used  in  scanning  to  MeV)  was  found  for  each 
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neutron  detector  and  for  each  day.  The  neutron  events  in  the 
chosen  time  interval  were  converted  to  an  MeV  scale  as  were 
the  accidental  neutron  events  and  the  neutron  events  accumulated 
during  positive  beam  running.  For  purposes  of  analysis  the 
events  were  grouped  into  intervals  of  0.2  MeV  (as  measured  on 
a  gamma  ray  scale) .  Several  corrections  were  applied  to  the 
neutron  events. 

a)  Accidentals;  Random  background  events  were  accumulated 
by  triggering  on  neutrons  occurring  during  the  15  p,sec  delayed 
time  gate.  For  the  3  and  4  in.  geometries  the  accidental  events 
rate  was  0.39  and  0.49  counts  per  p,sec  per  million  muon  stoppings, 
respectively.  As  has  been  pointed  out,  the  accidental  rate  in 
the  real  gating  interval  was  reduced  by  the  requirement  that  no 
anticoincidence  counter  pulse  (A-pulse)  appear  on  the  trace. 

The  effect  of  this  requirement  was  measured  by  considering  the 
time  distribution  of  A-pulses  relative  to  neutron  pulses  and 
also  with  respect  to  muon  pulses.  A  fraction  of  the  A-pulses 
were  correlated  in  time  with  neutron  pulses,  these  represent 
NA  failures;  others  appeared  earlier  in  time  than  the  muon 
pulse,  these  are  random  and  a  correction  is  to  be  applied  for 
having  discarded  these  few  events;  the  rest,  namely,  those  that 
appear  after  a  muon  but  uncorrelated  with  a  neutron,  are  the 
ones  of  interest  here  and  represent  reduction  of  accidentals. 

The  integral  of  the  energy  spectrum  of  accidental  neutron 
events  (Fig.  17)  was  corrected  for  the  length  of  the  "real" 
time  interval  and  for  the  above  reduction, and  the  normalized 
energy  spectrum  was  subtracted  point  by  point  from  that  of  the 


real  neutron  events . 
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b)  Neutron  Events  in  the  Positive  Beam;  in  order  to 
evaluate  possible  spurious  events  due  to  the  hig*i  flux  of  decay 
electrons  entering  the  neutron  detectors ,  it  was  desirable  to 
count  decay  electrons  in  the  absence  of  muon  capcure  neutrons. 
This  situation  is  readily  obtainable  by  stopping  positive  muons 
in  the  hydrogen  target.  Positive  muons  have  essentially  no 
probability  for  nuclear  interactions,  but  their  decay  electron 
energy  spectrum  and  angular  distribution  are  identical  with 
those  of  negative  muons.  It  was  found  that  a  non-negligible 
rate  of  neutron  events  persisted  under  positive  beam  running 
conditions.  These  events  were  found  to  have  the  time  distribu¬ 
tion  of  decay  electrons  and  an  explanation  for  t'  em  was  found 

in  photonuclear  processes  occurring  in  the  neutron  detectors. 

12  11  12  11 
The  photonuclear  reactions  C  (y,  n)  C  a:\d  C  (<y,p)  B 

show  large  increases  in  cross  section  between  20  and  30  MeV, 

29  30  31 

the  region  of  the  giant  dipole  resonance.  '  '  Somewhat 

less  probable  are  similar  reactions  involving  th**ee  particle 
final  states.  The  photoproton  reaction  has  a  higher  but  nar¬ 
rower  peak  in  its  cross  section  than  does  the  photoneutron 
reaction  which  reaches  a  maximum  at  22  MeV  and  falls  off  slowly. 
Between  zero  and  50  MeV  the  average  cross  section  for  the  two 
processes  are  approximately  equal.  The  energy  distribution  of 
the  proton  peaks  at  5.5  MeV  and  has  a  width  at  half  maximum  of 
2.8  MeV.  While  it  is  difficult  to  estimate  the  number  and 
energy  spectrum  of  gamma  rays  from  electron  bremsstrahlung  that 
are  available,  the  observed  rate  of  neutron  events  in  the  posi¬ 
tive  beam  is  roughly  consistent  with  photonuclear  processes 
occurring  in  the  thick  neutron  detectors. 
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The  possibility  of  stopping  positive  muons  in  the  hydrogen, 
where  the  physical  situation  is  identical  except  for  the  ab¬ 
sence  of  capture  neutrons ,  allows  an  accurate  correction  to  the 
data  independent  of  knowledge  of  the  source  of  events.  A  further 
confirmation  of  the  legitimacy  of  the  method  of  correction  is 
obtained  by  considering  the  neutron  events  in  the  energy  spec¬ 
trum  which  lie  above  the  muon  capture  cut-off  energy.  Such 
events  in  the  spectrum  obtained  from  the  negative  beam  are 
exactly  subtracted  by  similar  events  in  the  positive  beam  data? 
it  follows  that  the  subtraction  at  lower  energies  must  also  be 
of  the  correct  size. 

That  the  positive  beam  neutron  events  are  electron  asso¬ 
ciated  is  further  confirmed  by  the  time  distribution  of  these 
events  which  can  be  seen  in  Fig.  18  to  be  a  2.2  |j,sec  exponential. 
The  neutron  events  were  subdivided  into  the  usual  0.2  MeV  inter¬ 
vals  and  were  corrected  for  positive  beam  accidentals;  they 
were  then  normalized  to  take  into  account  the  number  of  decay 
electrons  counted  during  negative  beam  running  and  to  take  into 
consideration  differences  in  the  gating  time  and  were  subtracted 
point  by  point  from  the  neutron  events  accumulated  during  the 
negative  beam  experiment.  The  approximate  rate  of  positive 
beam  events  is  one  event  per  10,000  detected  decay  electrons, 
where  an  event  is  defined  as  a  proton  recoil  of  between  1  and 
9  MeV.  The  rate  within  the  energy  interval  of  interest  is 


of  course  lower. 
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c)  Wall  Capture  Events ;  The  neutron  time  distribution 
reveals  a  contribution  from  muon  capture  in  iron;  it  will  be 
shown  later  that  there  is  no  reason  to  expect  a  substantial 
contribution  from  any  material  other  than  iron.  The  energy 
spectrum  of  neutrons  from  muon  capture  in  iron  was  determined 
in  two  ways.  In  a  preliminary  measurement  a  study  was  made  of 
the  neutron  yields  and  spectra  from  muon  capture  in  various 
materials  including  carbon,  iron,  silver,  zirconium  and  lead. 
Figure  20  shows  the  spectra  obtained  from  these  materials  using 
a  neutron  detector  and  muon  beam  similar  to  the  ones  used  for 
the  principal  experiment.  An  additional  determination  of  the 
iron  spectrum  was  obtained  by  considering  only  the  early  time 
segment  of  the  actual  data  and  subtracting  from  this  the  appro¬ 
priate  number  of  accidentals  and  of  events  with  2.2  |j.sec  life¬ 
time.  The  energy  distribution  so  obtained  is  shown  in  Fig.  21; 
it  is  fitted  by  an  exponential  distribution  as  would  be  expected 
for  an  evaporation  source  but  is  even  more  closely  approximated 
by  an  exponential  at  low  energies  and  a  slowly  decreasing  func¬ 
tion  at  the  upper  end.  This  shape  is  expected  on  the  basis  of 
models  of  muon  capture  in  heavy  nuclei.  The  neutron  yield  from 
iron  was  normalized  to  the  number  obtained  from  the  time  dis¬ 
tribution  and  the  gating  interval  location  and  was  subtracted 
from  the  neutron  events  spectrum.  The  gating  interval  was 
chosen  so  as  to  keep  the  subtraction  small;  iron  events  would 
soon  overwhelm  the  data  at  earlier  times  and  the  subtraction 
is  insufficiently  reliable. 
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d)  Scanning  and  Statistical  Corrections :  The  maximum 
film  advance  rate  in  the  oscilloscope  photography  system  was 
1  frame  per  second.  It  was  considered  undesirable  to  employ  a 
blocking  time  of  this  magnitude  for  the  oscilloscope  trigger 
since  this  would  make  it  impossible  to  count  directly  those 
events  which  fell  within  the  dead  time.  Instead,  during  film 
scanning  an  IBM  card  record  was  made  of  events  which  overlapped 
and  which  were  unmeasurable . 

The  number  of  neutron  events  were  corrected  for  the  number 
of  events  lost  because  of  confusing  overlap  of  pulses  on  the 
trace  and  because  of  superimposed  traces?  correction  was  also 
made  in  the  categories  of  accidentals  and  positive  beam  events. 
The  percentage  correction  that  was  applied,  together  with  its 
uncertainty,  is  tabulated  below  for  each  type  of  loss. 


Events  with  2  mu  pulses 
Events  with  2  neutron  pulses 
Superposed  frames 
Accidental  pion  pulse 
Accidental  electron  pulse 


3  In. 

2.1%  ±  0.2% 
6.2%  ±  0.3% 
2.2%  ±  0.2% 
2.0%  ±  0.4% 
0.5%  ±  0.1% 


4  In. 

4.5%  ±  0.3% 
9.6%  ±0.3% 
4.8%  ±  0.3% 
1.6%  ±  0.3% 
0.8%  ±  0.1% 


Total  Correction 


13.0%  ±  0.6% 


21.3%  ±  0.6% 


This  correction  contributes  only  a  negligible  error  to  the 
final  result  by  virtue  of  a  complete  tabulation  and  record  of 


lost  events . 
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The  scanning  efficiency  demands  no  correction  since  essen¬ 
tially  every  film  frame  was  accounted  for  on  an  IBM  card,  and 
rejection  of  events  took  place  during  analysis  of  the  cards. 

Electronic  dead  times  are  all  too  short  to  have  an  appre¬ 
ciable  effect  on  the  rate  of  "events"  triggers  and,  since 
detector  pulses  are  displayed  directly,  can  have  no  effect  on 
the  pulses  appearing  on  the  traces. 

The  neutron  pulse  height  spectrum  was  plotted  for  each 
detector  in  each  geometrical  arrangement.  Detector  was  elimi¬ 
nated  from  detailed  consideration  early  in  the  analysis  since 
the  quality  of  neutron-gamma  ray  discrimination  was  doubtful 
for  a  large  fraction  of  the  data.  The  calculated  detector 
response  was  fitted  to  each  energy  spectrum  by  equalizing  the 
areas  above  some  chosen  energy  (on  the  gamma  ray  scale) ?  no 
scale  change  or  other  adjustment  was  needed  beyond  the  original 
conversion  of  arbitrary  pulse  height  units  to  energy  (MeV) 

using  source  data.  Figure  19  shows  the  six  plots  along  with  the 
2 

X  calculated  after  the  expected  curve  was  fitted. 

For  purposes  of  presentation  and  to  calculate  the  capture 
rate  the  neutron  data,  as  well  as  the  corrections,  were  grouped 
into  two  energy  blocks,  one  beginning  at  a  low  energy  depending 
on  the  detector  and  ending  at  2.8  MeV  (corresponding  to  5.8  MeV 
proton  recoils) ,  and  the  other  entirely  above  the  expected 
neutron  energy  from  3.2  to  4.2  MeV  (gamma  scale) .  The  data  in 
the  first  block  was  used  to  obtain  the  capture  rate  as  in 
Table  V,  while  the  second  block  in  Table  VI  demonstrates  the 
cancellation  of  events  above  the  expected  energy. 
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5.  Small  Modifications  to  the  Capture  Rate 

a)  The  presence  of  deuterium  in  the  hydrogen:  The  fraction 
of  muons  being  transferred  to  the  deuterium  branch  of  the  diagram 
in  Fig.  1  is 


cX 


F  = 


X  +  X  +  cX 
pp  o  e 


the  concentration  of  deuterium  in  the  target  is  less  than  one 
part  per  million.  Assume  the  measured  values  for  the  constants 
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then  F  =  0.007. 

Actually,  not  all  the  muons  which  transfer  to  deuterium  fail  to 
give  a  measured  neutron;  in  any  case  the  effect  of  deuterium  is 
negligible . 

b)  Other  impurities  in  the  hydrogen:  Assume  a  transfer 

rate  to  higher  atomic  number  impurities  equal  to  5  x  10^  sec 

-12 

the  concentration  of  these  impurities  is  ~  10  ,  certainly 

-9  -6 

less  than  10  .  In  this  case  F  =  25  x  10  ,  but  the  capture 

rate  is  higher  than  that  of  hydrogen  by  ~  200;  then  the  fraction 

of  neutrons  from  capture  in  these  impurities  is  0.005  at  most; 

again  negligible. 

c)  Captures  in  carbon  (outside  the  target) :  There  are 
only  three  possible  sources  of  carbon  capture;  great  efforts 
were  made  in  the  construction  of  the  target  to  eliminate  all 
low  Z  material,  and  plastic  scintillator  was  the  only  remaining 
potentially  dangerous  material. 
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i)  Muons  stopping  in  counter  #3a:  The  accidental  rate  of 
the  (123a) 3b  coincidence  is  10  per  minute.  These  accidentals 
could  very  likely  represent  stoppings  in  the  #3a  counter.  In¬ 
cluding  the  ratio  of  solid  angles  of  this  detector  to  the  hydrogen 
(0.06)  and  the  ratio  of  capture  rates  (100),  the  fractional  yield 

_3 

due  to  this  mechanism  is  2.5  x  10  ;  this  is  completely  negligible 

ii)  Muons  stopping  in  counter  #3b:  The  differential  stop¬ 
ping  distribution  curve  for  muons  (Fig.  9)  indicates  that  in 
the  vicinity  of  counter  #3b  the  stopping  rate  is  500  per  gram, 
while  in  the  peak  of  the  distribution  the  rate  is  3000  per  gram. 
The  counter  is  0.05  grams  thick  so  that  the  ratio  of  muons  stop¬ 
ping  in  the  counter  to  those  stopping  in  hydrogen  is  0.002. 

The  ratio  of  solid  angles  is  0.07,  the  extra  length  of  hydrogen 
causes  further  attenuation  of  the  neutron  yield  by  about  a  factor 
of  2,  so  the  percentage  neutron  yield  from  this  source  is  0.7%. 

iii)  Muons  penetrating  the  anticoincidence  detectors? 

This  case  is  an  important  one  to  study  carefully  since  the  solid 
angle  for  detection  of  any  neutrons  produced  in  the  anticoin¬ 
cidence  detectors  is  not  smaller  than  that  for  hydrogen  neutrons. 
It  will  be  possible  to  obtain  a  sufficiently  accurate  estimate 
of  this  source  of  background  by  considering  the  number  of  muons 
stopping  in  the  target  walls.  The  latter  number  can  be  obtained 
in  the  following  ways:  First,  the  iron  in  the  target  wall  com¬ 
prises  0.495  grams  while  the  higher  Z  components  together  com¬ 
prise  2.69  grams  (with  an  average  Z  of  42). 

Analysis  of  the  neutron  data  taken  on  film  gives  the  follow¬ 
ing  result  for  the  ratio  of  muons  stopping  in  iron  to  those 
stopping  in  hydrogen: 
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4  Inch  geometry  0.78%  ±0.11 

3  Inch  geometry  0.55%  ±  0.08. 

The  neutron  time  distributions  as  recorded  on  the  digitron 
clearly  show  both  a  high  Z  component  and  an  iron  component  in 
the  lifetime.  From  this  data  one  can  obtain,  independently  of 
solid  angle  and  efficiency,  the  ratio  of  muon  stoppings  in 
high  Z  to  the  stoppings  in  iron: 
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while  the  ratio  of  densities  is  5.45.  One  then  has  for  the 
percentage  of  muons  stopping  in  the  higher  Z  component  of  the 
walls : 

4  Inch  2.8%  ±  0.7 

3  Inch  3.4%  ±  0.8 

Second,  an  independent  estimate  of  the  number  of  muons 
stopping  in  high  Z  materials  comes  from  a  comparison  of  the 
yield  of  decay  electrons  per  muon  in  the  positive  and  negative 
beams;  from  this  source  one  obtains: 

4  Inch  5.5%  ±  2.2 

3  Inch  3.5%  ±  1.4 

This  is  in  rough  agreement  with  the  previous  estimate.  Analysis 

of  the  time  distribution  of  decay  electrons  at  early  times 

shows  an  apparent  11%  stopping  in  high  Z  component;  however, 

this  number  should  be  reduced  by  a  factor  of  2  or  3  because 

it  includes  neutrons  which  have  scattered  in  both  the  anti- 

coincidence  detector  and  the  neutron  detector  thus  simulating 

electrons;  this  is  important  in  the  case  of  high  Z  materials 

only  since  there  are  twenty  times  as  many  neutrons  as  electrons 
per  muon  stopping. 
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The  above  determinations  of  wall  stopping  yield  a  stopping 
rate  for  muons  in  the  vicinity  of  the  target  walls  of  1.4%  of 
the  muons  per  gram  of  wall  material.  Since  muons  are  at  the 
end  of  their  range  having  penetrated  much  hydrogen,  etc.,  the 
rate  of  stopping  is  surely  decreasing  with  increasing  thickness 
and,  therefore,  is  expected  to  be  somewhat  lower  than  1.4%  at 
the  position  of  the  anticoincidence  detector;  nevertheless,  this 
number  will  be  used  to  obtain  an  upper  limit.  A  further  check 
on  the  above  estimate  is  the  observation  that  the  electronic 
requirement  (123a3b)A  reduces  the  counted  number  of  muons  by 
(1  ±  0.3)%  below  the  (123a3b)  rate;  since  the  anticoincidence 
detector  is  1  gram  thick,  the  estimate  of  muons  penetrating  it 
is  verified. 

Analysis  of  film  taken  of  decay  electron  events  shows 
that  one  can  easily  measure  pulses  whose  size  is  one-fifth 
the  size  of  the  pulse  due  to  a  minimum  ionizing  particle  pass¬ 
ing  through  the  detector  (Fig.  22) .  For  the  anticoincidence 
detectors  of  the  thickness  used,  one-fifth  minimum  ionizing 
corresponds  to  an  energy  loss  of  0.5  MeV.  A  stopping  muon 
which  leaves  0.5  MeV  in  the  plastic  comes  to  rest  after  0.0035 
grams .  The  mylar  and  aluminum  coating  on  the  counters  are 
0.007  grams  thick.  Assuming  a  total  low  Z  thickness  of 
0.011  grams  wherein  a  stopping  muon  would  go  undetected,  one 

has  a  ratio  of  carbon  stopping  muons  to  muons  stopping  in 

-4 

hydrogen  1.5  x  10  .  Including  a  solid  angle  enhancement 

factor  and  a  capture  rate  factor,  we  have  as  an  upper  limit  on 
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the  neutron  yield  ratio  1.8%.  A  more  realistic  estimate  is 

probably  1%.  Although  it  is  much  smaller  than  the  quoted  error, 

it  is  included  as  a  correction  factor:  (-1  ±  0.5)%. 

d)  The  admixture  of  (p,p)  atomic  capture  in  the  data: 

The  (pp)  atomic  system  converts  to  the  (ppp)  molecular  ion  with 

6  —1 

a  rate  X  =  1.9  x  10  sec  or  a  partial  lifetime  of  about 

pp 

0.5  psec.  The  data  that  were  used  in  the  analysis  were  taken 
in  a  time  interval  beginning  at  1.2  psec  so  that  the  admixture 
of  (pp)  is  quite  small;  nevertheless,  it  is  easy  to  calculate 
the  exact  effect  on  the  capture  rate.  Using  the  usual  coupling 
constants,  the  capture  rate  in  the  molecular  state  is  expected 
to  be  560  sec  ^  while  the  singlet  atomic  capture  rate  is  about 
640  sec  ^ .  The  rate  in  a  mixture  of  the  two  which  is  observed 
from  t  =  tQ  to  infinity  is: 


R  =  X 
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In  the  present  experiment  the  effect  is  very  small  but  in  a 
bubble  chamber  experiment  ( t  =  0)  the  rate  is  expected  to  be 
higher  by  about  4%.  Of  course,  such  a  correction  presupposes 
knowledge  of  the  actual  rate. 
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e)  Possible  molecular  ortho-para  conversion:  The  rate  of 

transition  of  the  ortho  state  to  the  para  is  expected  to  be 

12 

quite  negligible  on  the  scale  of  muon  lifetimes.  One  can 

obtain  an  upper  limit  to  the  decay  rate  from  the  neutron-muon 

time  interval  distribution  curves  (Fig.  16)  by  noting  that  the 

capture  rate  in  the  para  state  is  expected  to  be  about  220  sec 

as  opposed  to  560  sec  ^  in  the  ortho  molecular  state.  The 

experimental  upper  limit  to  the  ortho-para  conversion  rate  XQ^ 

6  —1 

is  then  0.05  x  10  sec  ,  corresponding  to  a  partial  lifetime 
of  20  p.sec.  In  the  unlikely  circumstance  that  this  were  in 
fact  the  conversion  rate  the  apparent  capture  rate  would  be 
reduced  to  about  525  sec  ^ . 

6 .  Summary 

The  final  estimate  for  the  capture  rate  is  464  ±  42  sec  1, 
a  small  correction  has  been  applied  for  the  effects  mentioned 
in  the  last  section.  It  can  be  seen  (Table  V)  that  the  error 
is  principally  due  to  three  sources:  i)  The  various  subtrac¬ 
tions  to  the  spectrum,  ii)  The  number  of  muons  stopping  in  the 
hydrogen,  and  iii)  The  detector  efficiency  calculation.  The 
statistical  error  in  the  raw  number  of  detected  neutrons  is 
rather  small  but  the  above  three  categories  contribute  roughly 
4%,  5%  and  5%,  respectively. 

Further  improvement  in  the  precision  of  such  an  experiment 
would  clearly  require  reduction  of  uncertainties  in  a  multitude 
of  areas  and  is  likely  to  be  very  difficult. 
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VI.  PION  CAPTURE:  A  CONFIRMING  MEASUREMENT 
A.  Introduction 

It  is  extremely  desirable  to  apply  the  target,  detector 
geometry,  and  efficiency  calculation  to  a  similar  experiment 
which  is  somewhat  easier  to  carry  out  and  whose  result  is  pre¬ 
dictable.  A  natural  choice  for  the  check  experiment  is  pion 
capture  which  proceeds  by  two  almost  equally  likely  modes. 

When  a  negative  pion  comes  to  rest  in  hydrogen  it  immediately 

•  32 

interacts  with  a  proton.  Recent  measurements  show  that  the 

time  interval  between  a  pion  velocity  of  0.006  C  and  nuclear 

•  .  .  -12 

interaction  is  about  2.3  x  10  sec  in  liquid  hydrogen.  No 
molecular  effects  take  place  in  the  case  of  pions .  The  inter¬ 
action  modes  that  occur  are 


(a) 

7T  +  p 

-*  n  +  7T1 

(b) 

TT~  +  p 

-*  n  +  7 

Mode  (a)  is  inhibited  by  lack  of  phase  space  and  by  the  s-state 
interaction  to  the  point  that  mode  (b)  can  compete  with  it. 

The  neutron  energy  in  mode  (b)  is  8.88  MeV  while  the  one  emitted 
in  mode  (a)  has  0.4  MeV  and  is  undetectable  in  this  experiment. 
The  Panofsky  ratio  which  has  been  measured  many  times  is  defined 
as  the  ratio  of  probabilities  a/b  =  P.  The  quantity  measurable 
in  this  experiment  is  b/(a  +  b)  =  R.  The  expected  value  of  R 
is  then  R  =  1/(1  +  P)  =  0.396  using  the  recent  value  for  P.^ 
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B.  Description  of  Experiment 

The  pion  capture  experiment  can  be  performed  with  no  change 
in  the  experimental  arrangement  other  than  that  needed  to  stop 
a  known  number  of  pions  in  the  target.  This  was  accomplished 
by  reverting  to  the  normal  unpurified  meson  beam  situation  and 
installing  the  time-of -flight  system.  The  output  of  the  time- 
to-height  converter  was  sent  to  a  pulse  height  window  which 
selected  pions  in  the  beam;  the  signal  from  the  window  output 
was  then  put  in  coincidence  with  the  usual  "muon  stopping" 
pulse,  providing  a  trigger  which  was  free  of  muons  and,  more 
important,  of  electrons.  Since  pion  capture  is  a  prompt  process 
the  muon  triggered  delayed  gate  was  shortened  in  length  and  its 
delay  was  decreased  so  as  to  be  in  prompt  coincidence  with  the 
pion  stopping  pulse  (Fig.  23).  Because  of  the  relatively  high 
"events"  rate,  it  was  possible  to  take  data  on  a  kicksorter  as 
well  as  on  film.  On  two  occasions  during  the  muon  capture  run 
the  pion  experiment  was  set  up;  the  "real  events"  trigger  gated 
a  pulse  height  analyzer  on  which  was  displayed  the  neutron 
detector  pulse  height  output.  Data  were  recorded  on  film  in 
exactly  the  same  way  as  during  the  muon  capture  experiment; 
during  scanning  a  prompt  pion  stopping  pulse  was  required  to 
identify  an  event. 

Neutron-gamma  ray  discrimination  was  checked  using  the 
pion  capture  neutrons  and  gamma  rays  themselves,  rather  than 
a  radioactive  source,  to  ensure  that  the  correct  energy  interval 
would  be  considered.  Energy  calibration  was  possible  using  the 
4.4  MeV  carbon  gamma  ray  from  a  plutonium-beryllium  source 
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*  24 

(a  +  Be  -»  n  +  c  -*  C  +  7)  and  the  Na  source  for  the  lower  end 

of  the  scale.  The  high  energy  point  was  fixed  using  the  6.14  MeV 

16  16 

gamma  ray  from  the  following  reaction:  n  +  0  -»  N  +  p, 

N16  -  e“  +  016*  +  v,  016*  -  016  +  y  (6.14  MeV) .  016  (water) 

is  irradiated  inside  the  cyclotron  shielding  wall;  since  the 
beta  decay  has  an  8  sec  half-life,  a  high-speed  pump  is  used 
to  send  the  water  out  to  the  detector  house  and  back. 

There  is  such  an  abundance  of  pion  capture  neutrons  that 
even  if  the  detector  energy  response  were  unknown,  it  would  be 
self-determining  at  the  8.9  MeV  point.  In  fact,  the  high  energy 
point  obtained  in  this  way  helps  to  confirm  the  earlier  inde¬ 
pendent  measurement  of  the  scintillator  energy  response  (see 
Appendix  I) .  Only  a  single  neutron  detector  was  used  for  the 
pion  capture  measurement  in  order  to  reduce  the  events  rate  to 
a  tolerable  level  and  to  simplify  the  pulse  height  normaliza¬ 
tion  problem  on  the  kicksorter  (Fig.  24) . 

C.  Data  and  Results 

To  calculate  the  fraction  of  pions  captured  in  the  one  ob¬ 
servable  mode  it  was  necessary  to  obtain  the  actual  number  of 
pions  stopping  in  hydrogen.  Since  the  (123a3b)  signal  was 
gated  on  pions  alone  (on  time-of-f light  criterion)  it  was  suffi¬ 
cient  to  apply  the  same  factor  (the  ratio  of  muons  stopping  in 
hydrogen  to  beam  counts)  as  was  used  in  the  muon  experiment. 

Only  one  correction  to  the  neutron  spectrum  data  was 
needed;  it  is  expected  that  the  same  number  of  pions  stop  in 
wall  materials  as  did  muons  in  the  principal  experiment.  How¬ 
ever,  wall  captures  here  are  coincident  in  time  with  pion  arrival, 
as  also  are  the  captures  in  hydrogen.  Thus  wall  captures  play  a 
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relatively  larger  role  in  pion  capture  where  the  device  of 
delayed  counting  is  not  effective.  A  crude  estimate  of  the 
number  of  neutrons  expected  from  wall  interactions  was  made 
and  agrees  with  the  actual  correction  which  was  performed  in  a 
self-determining  manner.  In  the  neutron  energy  spectrum  plot 
a  flat  low  level  of  counts  persisted  beyond  the  clearly  visi¬ 
ble  cut-off  energy  of  hydrogen  events.  The  low  background 
level  was  extrapolated  back  toward  zero  with  the  appropriate 
exponential-like  rise.  The  total  subtraction  in  the  energy 
region  of  interest  amounted  to  12.8  1  3%.  The  neutron  energy 
distribution  obtained  from  analysis  of  the  film  was  plotted 
and  fitted  (by  normalization  only)  with  the  calculated  detector 
response  curve  (Fig.  25) .  Only  a  change  in  one  parameter  (the 
initial  neutron  energy)  was  needed  to  convert  the  Monte  Carlo 
calculation  of  detector  response  from  the  muon  capture  to  the 
pion  capture  situation.  Table  VII  lists  the  data  and  results 
of  three  independent  runs  of  the  pion  capture  experiment. 

TABLE  VII 

Calculation  of  Pion  Capture  Branching  Fraction 


Run 

Detected 

Neutrons 

Corrected 

Neutrons 

Pion  , 
Stops ( 10°) 

Actual  3 

Neutrons ( 10 j 

Branching 
Fraction  R 

a 

720 

628±6 . 4% 

0.19215. 5% 

77.818.2% 

0.405 

b 

1468 

128015. 4% 

0.387 

159.017.4% 

0.411 

c 

685 

597±7 . 3% 

0.185 

74.118.9% 

0.401 

Average  0.407±0.035 

Expected  0.396 
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The  error  of  the  mean  is  8.5%.  The  weighted  average  branching 
fraction  is  0.407  ±  0.035.  The  branching  fraction  implied  by 
the  accepted  value  of  the  Panofsky  ratio  of  1.53  is  0.396.  It 
should  be  pointed  out  that  while  the  statistical  errors  are 
independent  for  the  three  runs,  the  errors  assigned  to  the 
number  of  pion  stoppings  and  to  the  efficiency  calculation  are 
not.  This  was  taken  into  account  in  the  calculation  of  the  error 
of  the  mean. 

The  good  agreement  between  the  measured  value  of  the 
branching  fraction  and  the  one  expected  on  the  basis  of  more 
precise  experiments  lends  credibility  to  the  assertion  that 
in  the  muon  capture  experiment  the  muon  stopping  number  and 
wall  effect  were  assessed  correctly,  and  that  the  detector 
efficiency  calculation  represented  the  physical  situation  to 
within  the  error  assigned. 

VII.  RESULTS  AND  CONCLUSION 

The  present  muon  capture  experiment  yields  a  single  experi¬ 
mentally  determined  number;  the  capture  rate  in  the  (ppp)  molecu¬ 
lar.  If  it  were  not  for  iron  capture  neutrons  appearing  at 
early  times  it  would  have  been  possible  to  extract  a  meaning¬ 
ful  value  for  the  (pp)  atomic  capture  rate.  An  experiment  free 
of  high  Z  capture  could,  in  principle,  extract  the  triplet 
capture  rate  independently  of  efficiency  and  some  backgrounds 
by  observing  the  time  variation  of  the  capture  rate  at  early 
times  (assuming  the  singlet  rate). 
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The  corrected  experimental  result  for  the  molecular  cap¬ 
ture  rate  is 

♦ 

464  ±  42  sec  ^ . 

In  the  basis  of  the  best  estimates  for  the  coupling  constants 
the  expected  molecular  capture  rate  is 
562  ±  60  sec 

where  the  error  reflects  the  uncertainty  in  the  molecular  wave 
function  as  well  as  in  the  calculation  of  the  capture  rate  (see 
p.  13)  . 

There  are  several  possible  interpretations  of  a  discrepancy 
between  the  observed  and  predicted  rates. 

i)  The  present  value  of  the  molecular  wave  function  over¬ 
lap  could  be  in  error  by  as  much  as  10%.  A  calculation  with 

21 

improved  wave  functions  is  now  under  way. 

ii)  The  equality  of  the  electron  and  muon  axial  vector 

coupling  constant  is  implied  by  the  observation  of  the  ratio 

of  the  reactions  v  —  e  +  v  and  T  -  [i  +  v  •  At  the  momentum 

transfer  of  muon  capture,  the  axial  vector  coupling  constant  is 

2  2 

modified  by  a  term  of  order  (m^  /M^  )  which  has  a  negligible 
effect  on  gA-  Anomalous  behavior  of  the  momentum  dependence 
is  possible  but  is  considered  unlikely. 

iii)  The  conserved  vector  current  theory  identifies  the 
momentum  dependence  of  the  vector  coupling  constant  with  that 
of  the  nucleon  form  factor.  This  identification  allows  the 
coupling  constant  gv^  to  be  evaluated  at  the  muon  capture  mo¬ 
mentum  transfer.  Effects  of  the  conserved  vector  current  have 

12  12 

recently  been  clearly  observed  in  the  beta  decay  of  the  B  -N 
multiplet.^  Experiments  on  the  decay  tt+  -»  tt°  +  e+  +  v  are  in 
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progress,  and  some  have  already  measured  a  rate  in  rough  agree¬ 
ment  with  the  prediction  of  the  conserved  vector  current  theory? 
In  view  of  the  success  of  this  hypothesis  it  is  unlikely  that 
either  the  vector  coupling  constant  or  the  so-called  weak  magne¬ 
tism  terms  (Eqs.  15  and  16,  Sec.  II)  differ  from  the  predicted 
values . 

iv)  It  is  conceivable  that  time  reversal  invariance  fails 

36 

in  muon  capture.  If  this  is  the  case,  then  the  vector  and 
axial  vector  coupling  constants  could  have  a  relative  phase 
factor  (see  Sec.  II)  and  the  capture  rate  in  the  singlet  state 
is  reduced  below  what  would  be  expected  if  gv  =  -gA-  In  the 
vicinity  of  180°  the  phase  factor  is  quite  insensitive  to  the 
measured  capture  rate  since  the  rate  is  at  a  maximum  there. 

v)  The  calculation  of  the  induced  pseudoscalar  coupling 

g  could  be  in  error.  The  usual  result  is  g  /g  =  8,  but  there 
P  P  A 

are  several  indications  that  this  is  in  disagreement  with  obser¬ 
vation.  The  large  angular  asymmetry  of  neutrons  observed  in 

the  capture  of  polarized  muons  by  spin  zero  nuclei  is  consistent 

3  7  5c 

with  a  larger  g  /g  ratio.  The  recent  observation  of  muon 

p  A 

16 

capture  induced  transitions  in  0  is  also  consistent  with  a 
larger  ratio  but  the  interpretation  depends  on  nuclear  wave 
functions . 

Using  the  expression  for  the  muon  capture  rate  (Eq.  23, 

Sec.  II)  with  g  /g.  and  g./g  as  independent  variables  the 

p  A  A  V 

molecular  capture  rate  was  computed.  A  band  was  plotted  (see 
Fig.  26)  in  the  plane  of  these  variables  that  includes  those 
points  consistent  with  the  results  of  the  present  experiment. 
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This  experiment  can  be  interpreted  as  providing  a  direct 
measurement  of  gA/<Jv  assuming  that  gv  is  given  by  beta  decay 
experiments  and  the  Conserved  Vector  Current  theory.  It  can 
be  seen  that  if  the  expected  pseudoscalar  term  is  retained  then 
gA/gv  =  -1.05  ±  0.08.  But  if  the  expected  ratio  of  A  to  V  is 
used,  then  a  ratio  of  g^/g^  =  15.5  ±  3.5  is  required  in  order 
to  fit  the  data.  In  a  hydrogen  experiment  the  triplet  capture 
rate  is  far  more  sensitive  to  the  pseudoscalar  term  than  is 
the  molecular  rate  which  is  largely  singlet. 

It  is  clear  that  the  experiments  on  muon  capture  that 
have  been  performed  up  to  this  time  (see  Table  VIII)  do  not 
yield  conclusive  evidence  for  a  universal  v-A  theory. 

TABLE  VIII 

Muon  Capture  Results 

Result  Expected 


Experiment 

(sec  ■'") 

( sec 

Columbia  a 

515  ±  85 

562 

Chicago  13 

435  ±  100 

580 

CERN  C 

420  ±  75 

580 

This  Experiment 

464  ±  42 

562 

a  See  Ref.  7. 

v 

See  Ref.  22. 
c  See  Ref.  23. 
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APPENDIX  I 

Proton  Response  in  NE  213 
A.  Introduction 

When  minimum  ionizing  particles ,  that  is  particles  for 
ciE 

which  is  independent  of  energy,  traverse  a  scintillating 
material  the  light  output  per  centimeter,  ^  ,  is  constant. 

One  might  expect  that  the  generalization  of  the  above  observa¬ 
tion  to  non-relativistic  particles  would  be  the  following: 

dS  _  dE 
dr  A  dr  ‘ 

39  . 

This  is  not  the  case  however;  the  correct  expression  is 
given  by 

dS  a(dE/dr) 

dr  l+b(dE/dr) 

An  oversimplified  picture  for  the  extra  "saturation"  term  in 
the  denominator  requires  the  heavily  ionizing  particle  to 
create  around  its  path  in  matter  a  very  densely  ionized  "column" 
in  which  recombination  with  subsequent  light  output  is  quenched 
by  the  strong  local  electric  fields.  The  more  heavily  ionizing 
the  particle  the  more  severely  quenched  is  the  initial  (~  2  nsec) 
light  output  from  the  vicinity  of  the  column.  Recombination 
and  de-excitation  at  later  times  (~  100  nsec) ,  characteristic 
of  drift  velocities  and  column  diameter,  always  occur  and  are 
responsible  for  the  slowly  decaying  component  of  light  which 
remains  linear  with  energy  loss.  Particle  detection  systems 
employing  fast  phototubes  and  high  frequency  components  (dis¬ 
tributed  amplifiers,  fast  discriminators,  oscilloscopes,  etc.) 
are  likely  to  be  sensitive  primarily  to  the  fast  nonlinear 
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component  of  light  output.  Therefore,  the  use  of  organic 
scintillators  as  energy  spectrometers  for  heavy  particles 
(e.g.,  protons)  requires  that  the  non-linearity  of  response 
of  the  scintillator  be  accurately  known  for  the  particle  of 
interest.  It  has  recently  become  clear  that  the  scintillation 
efficiency  or  light  output,  at  least  in  inorganic  scintillators, 
depends  on  the  mass  of  the  particle  as  well  as  on  dE/dx.  This 
has  been  explained  by  considering  the  ionization  produced  by 
delta-rays  which  are  able  to  escape  from  the  primary  column. 

In  the  case  of  heavy  ions  a  substantial  fraction  of  the  total 
energy  available  goes  into  delta  rays.  Unfortunately  no  calcu¬ 
lations  are  available  at  this  time  for  organic  liquids . 

B.  Description  of  Experiment 
Experiments  were  carried  out  using  the  Nevis  Neutron 
Velocity  Selector  to  determine  the  response  of  the  organic 
liquid  scintillator,  NE  213,  to  neutron  induced  proton  recoils 
in  the  range  2  to  6  MeV.  In  this  energy  range,  in  a  hydrogen¬ 
ous  substance,  ionization  due  to  incident  neutrons  is  entirely 
from  elastic  hydrogen  collisions  (at  least  at  times  short  com¬ 
pared  to  thermalization  time) .  The  mean  free  path  for  4  MeV 
neutrons  in  NE  213  (density  =  0.88,  composition  =  CH.^  2)  is 
roughly  4  in.  Thus,  if  one  must  measure  the  proton  response 
function,  a  thin  detector,  in  which  multiple  collision  events 
are  improbable,  should  be  used;  otherwise  the  observed  response 
function  contains  an  implicit  dependence  on  detector  geometry 
and  initial  neutron  direction. 
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A  glass  cylinder  1/2  in.  thick  and  4  in.  in  diameter  was 
filled  with  oxygen-free  NE  213  and  oriented  coaxially  with  the 
neutron  beam;  it  was  viewed  by  a  5  in.  phototube  (EMI  9530  B) 
in  direct  contact  with  the  4  in.  face.  The  detector  was  placed 
at  the  end  of  a  helium  pipe  100  meters  from  the  source  of  a 
pulsed  neutron  beam.  The  beam  is  obtained  by  deflecting  the 
cyclotron  proton  beam  into  a  target  once  every  cyclotron  cycle 
(1/60  sec) .  The  uncertainty  in  the  initial  time  is  at  the  very 
worst  4  proton  rotation  times,  or  about  200  nsec;  at  best  it 
is  one  cycle,  or  50  nsec.  The  use  of  "time • of -flight "  tech¬ 
niques  permits  the  convenient  measurement  of  neutron  energy  and 
the  elimination  of  gamma  ray  background.  A  similar  experiment 
performed  near  the  neutron  target  of  a  Van  de  Graaf  generator, 
for  example,  would  require  efficient  neutron-gamma  ray  discrimi¬ 
nation  and  a  knowledge  of  the  discrimination  efficiency  as  a 
function  of  energy.  Pulses  directly  obtained  from  the  photo¬ 
tube  were  displayed  on  an  oscilloscope  and  photographed.  The 
oscilloscope  was  triggered  by  the  beam  deflection  pulse;  neutron 
time  of  arrival  was  measured  from  the  time  of  arrival  of  gamma 
rays.  The  pulse  height  output  was  calibrated  using  the  Compton 
recoils  from  gamma  ray  sources  of  known  energy,  and  checked  with 
cosmic  rays  traversing  the  detector  parallel  to  the  cylinder 
axis  (being  minimum  ionizing,  they  leave  a  calculable  amount 
of  energy  in  the  detector) . 

C.  Procedures 

1)  The  sinusoidal  output  of  a  5  Me  crystal  controlled 
oscillator  was  displayed  on  the  1  (j,sec/cm.  Tektronix  555  oscil¬ 
loscope  sweep  and  photographed  to  calibrate  the  sweep.  A  pulse 
train  consisting  of  one  fast  pulse  per  (xsec  was  photographed 
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to  check  that  the  vertical  deflection  sensitivity  was  not  a 
function  of  horizontal  position  on  the  sweep;  this  was  the  case 
to  better  than  3%.  These  photographs  were  scanned  under  exactly 
the  same  conditions  as  the  data,  and  thus  furnished  a  check  on 
the  reprojection  system  and  measuring  technique. 

2)  With  the  detector  in  precisely  the  same  position  and 

attitude  as  during  the  data-taking  run,  it  was  exposed  to  the 

24 

gamma  radiation  from  Na  source  which  was  prepared  by  irradi- 
23 

ating  Na  Cl  with  a  high  flux  of  neutrons  in  the  cyclotron. 

The  gamma  ray  energies  from  this  source  are  1.38  MeV  and  2.76  MeV. 

Several  hundred  Compton  recoils  were  photographed  and  scanned, 

and  a  pulse  height  spectrum  was  plotted.  The  interpretation 

of  the  spectrum  in  a  thin  detector  presents  some  difficulty 

since  total  absorption  of  energy  is  very  unlikely.  In  fact, 

a  single  Compton  collision  is  most  probable  and  the  upper  end 

of  the  spectrum  is  the  "Compton  edge";  the  secondary  gamma  ray 

is  left  with  at  least  0.2  MeV  which  might  not  be  converted  in 

the  detector.  As  a  check  on  the  procedure  for  obtaining  gamma 

ray  energy  points  the  pulse  height  spectrum  from  cosmic  rays, 

which  were  selected  by  a  coincidence  requirement  to  be  directed 

within  10°  of  the  vertical,  was  displayed  in  a  pulse  height 
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analyzer  and  compared  with  the  spectrum  from  Na  ;  the  energy 
loss  of  the  cosmic  rays  was  calculated  using  the  density  of  the 
scintillator,  etc.,  and  the  point  was  located  with  respect  to 
the  spectrum  of  the  source,  providing  an  extra  fiducial  energy. 

3)  The  tenth  dynode  of  the  photomultiplier  was  matched 
to  125  n  cable  through  an  emitter  follower,  the  signal  was  led 
directly  to  the  input  preamplifier  of  the  oscilloscope.  No 


-72- 


intermediate  circuitry  was  present  to  introduce  non-linearities. 

The  possibility  of  space  charge  saturation  in  the  photomulti- 
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plier  itself  was  avoided  by  studying  the  Na  spectrum  as  a 
function  of  photomultiplier  high  voltage  and  observing  that 
the  relationship  between  the  two  "edges"  of  the  spectrum  as 
well  as  the  shape  at  the  upper  end  remained  constant  for  at 
least  150  volts  above  the  1200  volts  chosen  as  an  operating 
point.  These  saturation  tests  were  performed  using  a  400  channel 
RIDL  pulse  height  analyzer  of  proven  linearity.  The  gamma  ray 
energies  of  the  source  encompass  the  neutron  energy  region  of 
interest . 

4)  Data  consisted  of  about  4000  photographs  of  the  oscil¬ 
loscope  sweep  showing  about  6  |isec  of  time  after  gamma  arrivals. 
On  a  second  sweep  the  beam  deflection  trigger  pulse  was  dis¬ 
played  but  was  used  only  as  a  trigger  and  as  a  crude  reference 
point  for  the  "zero-time"  gamma  pulse  with  respect  to  which 
neutron  times  were  measured.  Pictures  were  taken  at  the  rate 
of  1  per  second,  a  limitation  of  the  camera  film  advance  system. 
Source  pictures  were  also  taken  at  1  per  second  selected  at 
random  from  a  much  higher  source  rate.  Sufficient  D.C.  current 
was  available  at  the  photomultiplier  dynodes  to  ensure  that 
there  were  no  rate  dependent  effects  on  the  pulse  height. 

The  fundamental  equation  of  the  time-of-f light  technique 
is  dE/E  =  -2  dt/t  so  that  is  one  wants  10%  energy  resolution 
at  5  MeV  (t  =  3.2  jj.sec)  the  time  has  to  be  measured  to  5%  or 
0.16  p.sec.  By  measuring  pulse  times  on  film,  at  the  sweep 
speeds  used,  dt  is  about  0.0  5  y,sec,  so  that  the  limitation  in 
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neutron  energy  resolution  is  due  to  initial  time  jitter  (0.2  p,s) 
and  the  necessity  of  grouping  data  to  improve  statistics  (from 
0.35  MeV  to  0.7  MeV,  depending  on  energy). 

D.  Data  Analysis 

The  scanning  of  the  film  yielded  a  time  and  pulse  height 
coordinate  for  each  neutron  pulse  recorded.  The  useful  energy 
range  was  broken  up  into  eight  sub-intervals,  and  the  pulse 
height  spectrum  in  each  of  these  was  plotted. 

The  spectrum  of  proton  recoils  from  a  single  hydrogen 
collision  is  1/Eq  dE  for  E  <  Eq.  Since  the  light  output  is 
nonlinear  with  proton  energy,  there  is  a  rise  in  the  spectrum 
toward  low  energies .  A  very  good  approximation  to  the  proton 
response  function  is 

R(E)  =  AE  +  BE2 

at  least  below  about  6  MeV.  In  the  above  equation  E  is  the 
proton  energy  in  MeV;  the  analagous  response  to  electrons 
would  be 

Rg(E)  =  E 

This  means  that  a  proton  of  energy  E  gives  as  much  light  as  an 

2 

electron  of  energy  (AE  +  BE  ) .  On  the  basis  of  this  type  of 
response,  one  can  calculate  the  expected  proton  recoil  spec¬ 
trum  shape  in  the  approximation  of  a  single  collision.  The 
universal  spectrum  for  arbitrary  initial  neutron  energy  Eq  is 

^  (A2  +  4BE) =  N(E) 

Eo 

where  the  variable  E  is  now  the  linear  pulse  height,  or  gamma 
ray  energy.  This  function  is  plotted  in  Fig.  27a. 
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In  conjunction  with  a  Monte  Carlo  calculation  of  detector 
efficiencies  the  problem  of  a  thin  (1/2  in.  thick)  counter 
exposed  to  a  parallel  beam  of  5  MeV  neutrons  was  run.  As  expec¬ 
ted  the  results  indicated  about  10%  multiple  collisions;  the 
calculated  and  expected  spectra  are  shown  in  Fig.  27.  Reasonable 
assumptions  were  made  about  the  response  function;  detector 
resolution  was  included  in  the  next  phase  of  the  calculation. 

Thus  in  principle  one  can  determine  the  response  function  by 
fitting  the  experimental  points  corresponding  to  a  fixed  neutron 
energy  and  determining  the  parameters  A  and  B.  Data  of  high 
statistical  quality  would  be  needed  to  obtain  the  parameters 
in  this  way.  What  is  done  in  practice  is  to  determine  the  end 
points  of  the  spectra  for  each  of  the  eight  energy  blocks 
separately,  and  to  plot  these  against  the  incident  neutron 
energies;  this  gives  directly  the  parabolic  response  curve. 

This  method  is  valid  because  above  2  or  3  MeV  the  recoil  spec¬ 
tra  are  almost  flat  so  that  when  resolution  is  folded  in,  as 
it  already  is  in  the  experimental  data,  the  position  of  the 
point  at  half  height  determines  the  end  point  very  well.  If 
more  data  were  available  a  combination  of  both  methods  would 
increase  the  reliability  of  the  low  energy  points.  In  that 
case  one  would  also  first  "unfold11  the  resolution,  which  was 
hard  to  measure  accurately  in  the  detector  used.  Its  excellent 
light  collection  and  small  volume  made  the  resolution  as 
measured  with  gamma  rays  different  from  that  characteristic 
of  proton  recoils . 
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Figure  28  shows  the  response  function  obtained  in  this 

.  25 

experiment  as  well  as  that  used  by  Brooks.  These  two  agree 

within  the  errors  of  measurement.  They  are  both  in  disagree- 
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ment  with  response  curve  for  NE  213  recently  published;  that 
measurement  used  a  cylindrical  detector  with  the  neutrons  inci¬ 
dent  perpendicularly  to  the  axis  and  viewed  with  a  6810A  photo¬ 
multiplier.  The  rise  time  of  that  tube  is  about  2  nsec  while 
that  of  the  9530B  is  about  25  nsec  but  that  difference  should 
have  little  effect  on  the  pulse  height  measured.  That  measure¬ 
ment  is  in  the  direction  of  decreased  light  output  or  increased 
saturation.  It  is  known  that  the  coefficient  b  in  the  equation 

ff  =  a(dE/dr)  (1  +  b  If)'1 

increases  as  the  amount  of  Terphenyl  in  the  scintillator  de¬ 
creases.  It  may  be  that  differences  in  composition  account  for 
the  discrepancy. 

Additional  data  of  high  statistical  accuracy  were  obtained 
by  gating  the  input  to  a  pulse  height  analyzer  and  displaying 
the  proton  recoil  spectrum.  The  delay  of  the  narrow  time  gating 
pulse  determined  the  neutron  energy  that  was  being  analyzed; 
its  width  determined  the  energy  resolution.  A  series  of  energy 
points  up  to  about  6  MeV  were  analyzed;  at  each  energy  setting 
gamma  ray  sources  were  used  to  check  the  linearity  of  the  gating 
system.  An  oscilloscope  was  used  to  monitor  the  time  interval 
between  the  arrival  of  the  gamma  ray  burst  from  the  cyclotron 
and  the  position  of  the  gating  pulse  which  selects  neutrons 
assuring  accurate  energy  measurement.  The  response  points 
obtained  in  this  way  are  shown  on  the  plot  (Fig.  28)  along  with 
those  obtained  from  analysis  of  the  film. 
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E.  Conclusion 

The  proton  response  curve  obtained  in  this  experiment 

25 

agrees  rather  closely  with  that  of  Brooks  .  A  simple  analytic 
fit  to  the  curve  is  obtained  by  assuming  a  parabolic  shape 
below  6  MeV  and  linear  above  that.  The  functions  that  were 
used  were : 

R  =  0.28  E  +  0.033  E2  E  <  6.  MeV 

R  =  0.70  E  -  1.330  E  >  6.  MeV 

The  response  curve  determines  the  shape  of  the  recoil 
spectrum  in  two  principal  ways ;  mainly  it  determines  the  upper 
energy  fall-off  point,  that  is,  the  response  to  a  neutron  which 
loses  all  its  energy  in  a  single  collision;  secondly  it  determines 
the  curvature  of  the  response  spectrum  near  the  low  energy  end. 
When  fitting  actual  data  the  first  characteristic  is  essentially 
self-calibrating,  this  means  that  even  a  small  error  in  the  end¬ 
point  would  result  in  a  large  deviation  of  the  curve  from  the 
data.  This  statement  is  weakened  somewhat  by  the  finite  resolu¬ 
tion  of  the  detectors,  nevertheless  a  5%  error  in  the  endpoint 
would  be  apparent.  The  second  characteristic  of  the  response 
curve  causing  the  spectrum  to  rise  near  zero  energy  is  difficult 
to  measure,  but  if  the  data  to  be  fitted  do  not  extend  down  into 
the  region  where  the  response  is  rising  sharply  then  the  inte¬ 
grated  yield  is  not  sensitive  to  a  precise  knowledge  of  the 
energy  scale  in  the  low  energy  region.  In  the  muon  capture 
experiment  the  lower  energy  cut-off  that  was  chosen  did  in  fact 
eliminate  the  steeply  rising  part  of  the  spectrum;  it  should 
be  noted  also  that  in  the  case  of  the  hydrogen  experiment  part 
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of  the  spectrum  increase  near  zero  is  due  to  neutron  degrada¬ 
tion  in  the  hydrogen  target,  so  that  the  actual  shape  is  some¬ 
what  less  sensitive  to  the  response  curve  than  normally.  A 
second  feature  of  the  hydrogen  experiment  that  affects  the 
rise  in  the  spectrum  near  zero  is  the  necessity  of  subtracting 
iron  events  and  accidental  events  from  the  neutron  spectrum; 
both  of  these  corrections  are  characterized  by  the  fact  that 
their  spectrum  increases  sharply  at  low  energies. 
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APPENDIX  II 

Measurement  of  Muon  Stopping  Distribution 

Both  the  solid  angle  subtended  by  the  neutron  detectors 
and  the  counting  efficiency  for  neutrons  (due  to  hydrogen 
self-absorption)  depend  on  the  spatial  distribution  of  muons 
that  have  come  to  rest  in  the  hydrogen  target.  Upon  the  com¬ 
pletion  of  the  muon  capture  experiment,  using  the  identical 
geometrical  arrangement  and  beam,  a  study  was  made  of  the  muon 
distribution  to  bring  up  to  date  preliminary  measurements  that 
had  been  made  earlier. 

A  dummy  target  consisting  of  a  detachable  replica  of  the 
front  target  window,  a  replica  of  the  side  walls  and  a  series 
of  alternating  layers  of  styrofoam  and  polystyrene  (which  simu¬ 
lated  the  stopping  power  of  liquid  hydrogen)  took  the  place  of 
the  actual  target.  The  pseudo-hydrogen  could  be  inserted  or 
removed  in  multiples  of  1/2  in.  thichnesses,  so  that  the  beam 
could  be  scanned  at  various  positions  along  the  target  axis. 
Scanning  was  accomplished  with  a  1"  x  1"  plastic  scintillator 
mounted  on  a  1"  diameter  RCA  (6199)  photomultiplier  which  was 
mounted  on  a  delicately  balanced  engine.  This  allowed  the 
detector  to  move  in  two  directions  perpendicular  to  the  beam 
as  well  as  along  the  beam  axis.  A  1/16  in.  aluminum  cap  was 
provided  to  fit  over  the  plastic  scintillator;  at  each  setting 
of  the  detector  position  the  beam  was  counted  with  the  cap 
alternately  on  and  off  to  obtain  the  number  of  muons  stopping 
per  centimeter  of  hydrogen. 
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Horizontal  and  vertical  scans  were  recorded  with  the  detec¬ 
tor  at  the  front  of  the  "target"  and  with  the  front  window  re¬ 
moved.  The  total  number  of  mesons  entering  the  target  region 
was  determined  by  computing 

2tt  f*rp(r)dr 

where  p(r)  is  the  average  radial  distribution  of  the  beam  ob¬ 
tained  from  the  two  scans.  This  number  was  compared  with  the 
number  of  muon  stopping  signals  as  defined  in  Sec.  V-B-l.  The 
ratio  of  actual  muons  to  "mu  stops"  in  the  4  in.  geometry  was 
found  to  be  0.75  ±  0.05  which  is  in  rough  agreement  with  the 
ratio  as  determined  from  decay  electrons  (see  Sec.  V-B-l) . 

The  error  due  to  the  finite  resolution  of  the  scanning  counter 
was  estimated  to  be  small  (2%)  compared  to  statistical  errors. 

Scans  were  performed  at  various  axial  positions  (Z)  in 
order  to  determine  the  radial  distribution  as  a  function  of 
depth  and  to  obtain  the  stopping  distribution  along  the  axis . 

In  order  to  improve  statistical  accuracy  it  was  necessary  to 
use  the  total  number  of  muons  [N(z)  =  2i7Trp(r)dr]  at  each  axial 
point;  the  distribution  obtained  in  this  way  is  then  equivalent 
to  the  one  obtained  by  range  measurement. 

The  radial  distributions  obtained  along  the  target  axis 
are  shown  in  Fig.  29.  No  differences  in  distribution  were  ob¬ 
served  in  the  3  in.  and  4  in.  geometrical  situations.  It  is 
known  that  the  number  of  wall  stopping  muons  is  different  in 
the  two  cases  (see  Sec.  V-B-5)  but  the  statistical  accuracy 
was  insufficient,  especially  in  the  wings  of  the  distribution, 
to  determine  that  number  from  the  scanning  data. 
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For  use  in  the  neutron  detector  efficiency  calculation 

several  simplifying  assumptions  about  the  distributions  were 

made.  The  spatial  stopping  distribution  was  taken  as  the 

product  F^(r)  •  F2(z)  and  both  functions  were  assumed  to  be 

gaussian.  The  standard  deviations  (o  and  az)  were  obtained 

from  the  experimental  curves  and  range  curves:  a  =  2"  , 

a  =  3.3"  . 
z 
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APPENDIX  III 

Calculation  of  Detector  Efficiency  and  Solid  Angle 
A  Monte  Carlo  calculation  designed  to  run  on  the  IBM  7090 
was  performed  in  order  to  calculate  the  neutron  detector  effi¬ 
ciency.  It  traced  neutron  trajectories  from  their  starting 
point  in  the  hydrogen  target  through  the  walls  and  anticoinci¬ 
dence  counter  and  into  the  neutron  detector  where  proton  recoils 
were  tabulated  to  yield  the  detector  response. 

The  geometrical  features  that  were  built  into  the  calcula¬ 
tion  included  the  cylindrical  hydrogen  target,  high  Z  target 
wall,  anti-coincidence  counter,  neutron  detector  walls,  neutron 
detector  volume,  lead  shielding  surrounding  detectors  and  gaps 
where  appropriate.  To  specify  the  geometry  and  to  make  the 
calculation  manageable  the  entire  domain  was  subdivided  into 
80  contiguous  convex  regions  bounded  by  natural  surfaces  when¬ 
ever  possible.  The  material  composition  of  each  region  was 
specified  by  the  atomic  weight  of  each  component  and  by  the 

atomic  density  of  each.  Neutron  elastic  and  inelastic  cross 

40 

sections  (supplied  by  NDA  )  were  tabulated  at  energy  intervals 
of  0.05  MeV  for  all  the  relevant  elements  and  were  used  to 
tabulate  neutron  mean  free  paths  for  each  type  of  material 
present.  Additional  parameters  were  tabulated  as  a  function 
of  energy  and  material;  these  specified  the  probabilities  of 
scattering  on  a  given  element  within  a  material  and  the  relative 
probability  of  inelastic  and  elastic  collisions.  It  was  assumed 
that  when  a  neutron  undergoes  an  inelastic  collision  it  is  lost. 
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The  distribution  of  neutron  starting  points  was  taken  to 

be  the  product  of  gauss ians  discussed  in  Appendix  II?  a  gaussian 

distribution  was  simulated  by  noting  that  the  average  of  a 

large  number  of  independent  randomly  distributed  numbers  on  a 

finite  interval  tends  to  a  normal  distribution.  The  starting 

neutrons  were  given  a  random  direction  in  space  and  a  specified 

initial  energy.  A  random  number  was  chosen  according  to  an 

exponential  distribution  by  a  technique  discussed  by  Von  Neu- 
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mann  ;  this  number  represented  the  number  of  mean  free  paths 
the  neutron  would  traverse  before  undergoing  a  collision.  The 
mean  free  path  in  the  region  was  then  compared  with  the  distance 
to  the  region  boundary  along  the  direction  of  particle  motion 
to  determine  whether  the  neutron  was  going  to  leave  the  region 
or  suffer  a  collision.  In  the  former  case  the  new  region  to  be 
entered  was  found  and  the  process  was  repeated;  in  the  later 
case  the  collision  type  was  decided  and  the  outgoing  energy 
and  direction  were  calculated.  If  the  neutron  energy  were 
above  a  certain  lower  cutoff  usually  taken  as  0.1  MeV,  then  the 
original  process  was  repeated;  if  not,  the  total  neutron  detec¬ 
tor  response  for  the  event  was  calculated  and  a  new  history  was 
begun.  When  a  hydrogen  collision  occurred  in  the  neutron  detector 
region  the  proton  recoil  energy  as  well  as  the  neutron  energy 
were  recorded;  at  the  end  of  a  history  when  the  neutron  either 
fell  below  the  cut-off  energy  or  left  the  domain  of  interest, 
the  responses  to  the  proton  recoils  of  that  history  were  summed 
and  the  resulting  response  was  tabulated.  Knowledge  of  the 
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proton  response  function  for  the  liquid  scintillator  used  was 
needed  to  combine  in  a  meaningful  way  the  individual  proton 
recoils  of  a  given  history.  After  a  determined  number  of  his¬ 
tories  was  completed  a  gaussian  energy  resolution  was  folded 
into  the  detector  response.  The  gaussian  width  was  taken  to 
vary  as  1//E  in  accordance  with  the  assumption  that  resolution 
broadening  is  due  to  statistical  fluctuations  in  photo-electron 
number.  Also  tabulated  were  histories  in  which  the  neutron 
produced  a  proton  recoil  in  the  anticoincidence  detector;  these 
were  subtracted  from  the  final  expected  yield  since  they  would 
be  electronically  rejected. 

It  was  possible  to  improve  the  statistical  accuracy  of  the 
Monte  Carlo  calculation  of  detector  response  by  a  device  called 
"splitting".  Since  only  about  20%  of  the  starting  neutrons 
were  destined  to  enter  the  neutron  detectors  the  remaining  80% 
of  the  histories  in  no  way  contribute  to  a  knowledge  of  the 
actual  response;  they  only  serve,  in  a  manner  of  speaking, 
to  specify  the  solid  angle  and  the  interactions  within  the 
hydrogen  target.  It  is  wasteful  to  spend  too  great  a  fraction 
of  the  calculation  time  on  those  lost  neutrons  so  when  a 
neutron  is  about  to  enter  the  detector  region  5  "latent  neu¬ 
trons"  (for  example)  are  produced  with  the  energy  and  spatial 
properties  of  the  parent  neutron.  These  latent  neutrons  are 
followed  in  turn  through  the  detector  and  any  response  that 
might  occur  is  weighted  by  one-fifth  when  it  is  added  to  the 
total;  thus  improved  statistics  are  obtained  on  the  random 
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processes  which  take  place  in  the  detector.  Neutrons  which 
have  entered  the  detector  and  which  then  leave  are  not  split 
if  they  happen  to  return  a  second  time. 

The  statistical  accuracy  of  the  calculation  was  estimated 
by  tabulating  results  at  the  end  of  each  of  a  series  of  short 
calculation  runs  (with  independent  random  numbers)  and  using 
the  fluctuation  as  a  measure  of  the  accuracy.  The  accuracy  of 
the  total  yield  is  of  course  greater  than  the  accuracy  of  the 
spectrum  shape.  The  main  calculation  was  performed  with  15,000 
starting  neutrons;  the  statistical  accuracy  of  the  integrated 
yield  was  about  2%,  which  is  smaller  than  the  estimated  error 
due  to  uncertainties  in  other  calculational  parameters. 

The  calculation  was  performed  with  several  initial  neutron 
energies  from  2.5  MeV  to  8.9  MeV  in  order  to  check  behavior  and 
for  possible  use  in  related  experiments.  The  relative  efficiency 
is  plotted  in  Fig.  30  as  a  function  of  neutron  starting  energy 
and  detector  energy  cut-off.  The  curve  shapes  are  typical  of 
those  of  similar  detectors  despite  the  presence  of  the  hydrogen 
and  the  shielding  materials . 

The  calculation  was  performed  under  a  variety  of  conditions 
to  teat  the  effect  of  various  features  of  the  problem.  Tests 
were  made  with  various  source  distributions  to  check  the  sensi¬ 
tivity  to  error  in  the  muon  stopping  distribution.  It  was 
found  that  if  the  width  of  the  axial  distribution  is  decreased 
by  about  12%,  then  the  yield  changes  by  about  3%;  the  yield  is 
about  half  as  sensitive  to  increases  in  axial  width.  For  a 
change  in  the  radial  width  of  about  20%  the  yield  changed  by 
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less  than  2%.  Uncertainty  in  measurement  of  the  geometrical 
arrangement  might  account  for  2%  or  3%  error;  fortunately  the 
symmetry  of  the  detector  locations  results  in  some  cancellation 
of  position  errors.  Geometrical  approximations  have  little 
effect. 

Although  the  presence  of  the  liquid  hydrogen  does  not 
change  the  integrated  neutron  yield  very  much  ( ahout  10%  reduc¬ 
tion)  it  has  an  important  effect  on  the  neutron  pulse  height 
spectrum.  Only  about  45%  of  the  neutrons  leave  the  hydrogen 
unscattered;  the  remainder  have  an  almost  flat  energy  spectrum 
from  zero  to  the  maximum  energy.  Thus  the  original  5.2  MeV 
monochromatic  neutron  is  badly  degraded  and  there  is  a  pile- 
up  of  events  with  small  pulse  heights  many  of  which  are  lost 
below  the  lower  energy  threshold.  The  8.9  MeV  neutrons  from 
pion  capture  are  degraded  to  a  lesser  extent  and  their  pulse 
height  spectrum  shows  the  characteristic  plateau. 

The  calculated  neutron  pulse  height  spectrum  was  folded 
with  several  resolution  functions;  if  the  energy  cut-off  point 
is  not  too  near  the  end  of  the  spectrum  different  resolution 
widths  have  a  negligible  effect  on  the  integrated  yield. 

The  calculation  that  has  been  described  was  carried  out  on 
an  IBM  7090  Digital  Computer.  The  computer  code  was  written 
in  the  form  of  a  logical  control  program  with  sixteen  smaller 
sub-programs  to  carry  out  specific  calculations  and  routines. 
Most  of  the  programs  were  written  in  the  Fortran  language  but 
several  frequently  used  segments  involving  random  number  calcu¬ 
lations  were  coded  in  the  computer  language  (FAP) .  The  geome- 
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trical  specifications  and  the  element  cross  section  tables 
were  stored  on  magnetic  tape;  parameters  such  as  initial  energy, 
starting  distribution  widths,  cut-off  energy,  total  number  of 
histories,  etc.  were  specified  each  time  the  calculation  was 
performed.  While  it  was  impossible  to  store  the  details  of 
each  history,  tabulations  of  intermediate  stages  in  the  neutron 
history  were  available  as  were  the  neutron  detector  response 
without  resolution  broadening. 

By  letting  the  materials  in  all  regions  except  the  neutron 
detectors  be  replaced  by  vacuum  one  could  approximate  the  solid 
angle  subtended  for  the  detection  of  decay  electrons;  in  addi¬ 
tion  the  tabulation  of  electron  path  lengths  in  traversing  the 
neutron  detectors  gave  an  estimate  for  the  edge  effects  in 
counting  electrons.  The  mean  depth  of  the  neutron  detector 
is  about  3.4  in.  for  particles  which  enter  the  detector  and 
which  originated  with  the  muon  stopping  distribution  in  the 
target;  the  most  probable  depth  is  about  4.6  in.  These  numbers 
indicate  that  the  mean  inclination  for  particles  which  enter 
the  detector  is  about  30  deg  from  the  perpendicular  to  the 
target  axis  and  that  the  effective  wall  thickness  is  about 
15%  larger  than  the  actual  thickness. 

The  calculated  solid  angle  for  the  four  neutron  detectors 
was  19.6%.  The  effective  solid  angle  for  electrons  was  taken 
to  be  18.5%;  the  reduction  was  made  to  exclude  electrons  which 
traversed  less  than  1/2  cm  of  scintillator  and,  therefore, 
deposited  less  than  1  MeV  of  energy  in  the  detector. 
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APPENDIX  IV 
Molecular  Spin  States 

The  (M-P)  atom  is  in  a  singlet  state  by  the  time  it  is 
likely  to  form  a  {p(j,p)  molecule.  It  has  been  shown  that  it 
is  most  probable  that  an  ortho  molecule,  with  spin  and  orbital 
angular  momentum  equal  to  1  is  formed  rather  than  the  ground 
state  para  molecule.  The  (pp,p)  system  then  has  the  following 
properties : 

Total  angular  momentum  (j)  =  1/2,  3/2 

Proton  spin  (S^)  =  1 

Total  spin  (S)  =  1/2 

It  has  been  shown  by  calculation  that  there  is  no  appreciable 
interaction  between  the  total  spin  and  the  orbital  angular 
momentum,  and  that  the  total  spin  is  1/2.  Under  these  condi¬ 
tions  one  can  calculate  the  relative  amounts  of  (pip)  singlet 
and  triplet. 

Consider  the  system  of  the  three  spins  s  ,  s  ,  and  s 

and  label  them  by  j^,  j ^ »  and  j^.  Define  j1  =  3 2  +  -*3  an<^ 

j1'  =  ,  then  the  ortho  molecular  state  is  characterized 

by  j 1  =1  and  it  is  desired  to  convert  to  the  representation 

where  j11  is  either  0  or  1.  A  unitary  transformation  connects 

42 

these  representations  : 


where  the  R's  are  related  to  the  Racah  coefficients  W  as  follows: 


Rjnj,  =  [  ( 2  j  "  +  1)  ( 2  j  '  +  l)]t  W(j1j2jj3;  j  •  j ") 
j  =  h  +  h  +  j3 
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The  coefficients  needed  are  10)  and  W(^^;  11)  ;  these 

are  tabulated  and  one  ends  with 

(R01)2  =  3/4  and  (Ru)2  =  1/4 

which  are  the  relative  weights  of  the  two  states .  In  the  case 

2 

of  the  para  molecule  ( Rq ^ , )  is  given  by  ( 2 j 1  +  l)/4.  Then 
( 2 j ’  +  1) /4  =  1/4  or  3/4  for  the  singlet  and  triplet  states, 
respectively. 

These  results  could  actually  have  been  obtained  with  more 
elementary  techniques,  but  Racah  coefficients  furnish  a  natural 
means  of  changing  representations  when  three  spins  are  coupled. 
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TABLE  II 

Target  Wall  Composition 


Material 

rv 

Side 

Wall 

Zt 

Thickness 
(in. ) 

Thickness 
( gm/cms ) 

Thickness  in 
Radiation  Length 

Silver 

47 

0.020 

0.535 

0.060 

Stainless 

Steel 

26 

0.025 

0.49  5 

0.033 

Silver 

47 

0.020 

0.535 

0.060 

Zirconium 

40 

0.100 

1.620 

0.155 

Total 

3.185 

0.308 

1 - - 

Front  Window 

Stainless 

Steel 

26 

0.013 

0.258 

0.017 

Silver 

47 

0.001 

0.027 

0.003 

Stainless 

Steel 

26 

0.010 

0.198 

0.013 

Total 

0.483 

0.033 
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TABLE  III 

Properties  of  Beam  Counters 


Photo- 

Size 

Thickness 

Number 

Multiplier 

( in. 

1_ 

(in.) 

Purpose 

0 

56AVP 

6  x 

6 

0.2500 

Time-of-flight  only 

1 

6810A 

8  x 

8 

0.0625 

Beam  telescope 

2 

56AVP 

6  x 

6 

0.3750 

Time-of-flight  and  ] 

3a 

6810A 

5  x 

5 

0.0625 

Beam  telescope 

3b 

6810A 

5  x 

5 

0.0200 

Final  Beam  Counter 

4 

6810A 

12  x 

12 

0.3750 

Anticoincidence  for 

Range  Study  Only 
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TABLE  V 

Calculation  of  Capture  Rate 
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TABLE  VI 

Calculation  of  Rate  of  Events  above  Muon  Capture 
Neutron  Energy 


Energy  Interval  -  3.2  -  4.2  MeV  (gamma  ray  scale) 


Neutron 

Counter 

4 

In.  Geometry 

3 

In.  Geometry 

N1 

N2 

Nl 

N2 

N4 

Neutrons 

17.0 

±  4.6 

12.1 

17.0 

18.1 

15.9 

13.6 

Accidentals 

2.9 

1.7 

1.2 

2.2 

2.8 

0.9 

0.0 

Positive  Beam 
Events 

23.8 

6.0 

5.9 

18.5 

14.6 

8.2 

14.8 

Iron  Captures 

1.0 

0.6 

0.8 

1.0 

0.7 

1.0 

Corrected 

Events 

10.7 

db  7.8 

4.4 

-4.5 

-0.3 

6.1 

-2.2 

Average 

-1.2 

±  3.0 
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FIGURE  CAPTIONS 

Fig.  1  -  A  diagram  showing  the  history  of  a  muon  after  coming 

to  rest  in  pure  liquid  hydrogen.  XQ  is  the  free  muon 

decay  rate,  Xc  is  the  muon  capture  rate  in  the  p.p 

atom  X1  _  is  the  molecular  capture  rate.  X  is  the 
cap  c  pp 

molecular  ion  formation  rate.  When  a  small  amount  of 
deuterium  is  present  (concentration  c)  then  the  process 
indicated  by  a  dashed  line  can  occur. 

Fig.  2  -  a)  The  four  fermion  interaction  responsible  for  muon 
capture  ignoring  strong  interactions. 

b)  An  intermediate  pion  state  diagram  identified  by 
the  conserved  vector  current  theory  as  due  to  a  basic 
pion-lepton  interaction. 

c)  The  contribution  to  the  "induced  pseudo-scalar" 
term.  Part  of  the  diagram  is  responsible  for  pion 
decay. 

d)  A  correction  term  to  the  axial  vector  current. 

Fig.  3  -  A  diagram  of  the  hydrogen  target  and  associated  puri¬ 
fication  system.  The  deuterium  is  used  for  the 
experiment  that  measures  molecular  constants. 

Fig.  4  -  The  hydrogen  target  and  some  local  shielding  shewing 
locations  of  detectors. 

Fig.  5  -  Cyclotron  floor  plan  showing  shielded  "house"  and 
beam  detectors. 

Fig.  6  -  Differential  beam  stopping  distributions  in  the  normal 
mode  of  operation  (114  MeV/c)  and  in  the  purified 
muon  beam  (132  MeV/c).  The  pion  range  is  between 
1  and  2  in.  of  carbon?  the  muon  range  is  3  to  4  in. 
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Fig.  7  -  Pulse  height  distribution  of  time  to  height  converter 

output.  Beam  particles  traverse  a  6  meter  flight  path. 
Lowest  curve  shows  the  unpurified  beam  with  electrons 
on  the  left,  muons,  then  pions.  Upper  curve  is  puri¬ 
fied  muon  situation. 

Fig.  8  -  a)  and  b)  Beam  path  before  hydrogen  target  showing 

beam  detectors,  collimators,  and  moderator.  The  "3  in. 
and  "4  in. "  geometries  were  used  in  successive  weeks 
of  running. 

Fig.  9  -  Stopping  distribution  of  purified  muon  beam,  showing 
number  stopping  per  gram  of  carbon  as  a  function  of 
moderator  thickness?  and  total  number  of  stopping 
particles . 

Fig.  10-  Counting  rate  in  decay  electron  telescope  as  beam 
stopping  point  is  shifted  within  hydrogen  target. 

Number  of  muons  in  beam  telescope  is  shown  also  as 
a  function  of  moderator  thickness  (in  inches  of  carbon) 

Fig.  11-  Circuit  used  for  neutron-gamma  ray  discrimination. 

Typical  pulse  shapes  are  shown;  pulses  due  to  gamma 
rays  are  shown  as  dashed  where  they  differ  from  those 
due  to  neutrons . 

Fig.  12-  a)  Photograph  of  live  oscilloscope  display  of  neutron 
or  gamma  ray  pulse  height  versus  slow  component  resi¬ 
due.  Horizontal  axis  is  energy,  increasing  toward  the 
right.  The  neutron  residue  is  displayed  on  the  verti¬ 
cal  axis  increasing  downward.  Upper  photograph  shows 
gamma  ray  and  neutron  "bands"  due  to  a  Pu-Be  source. 
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Lower  picture  is  "band"  due  to  pure  gamma  ray  source 

1^7 

(Cs  -  0.67  MeV)  . 

b)  Neutron  and  gamma  ray  "bands"  for  a  Pu-Be  source. 
Lower  picture  shows  effect  of  a  residue  discrimina¬ 
tor  used  to  select  only  neutrons. 

c)  and  d)  Plots  of  pulse  height  versus  residue  for 
neutron  and  gamma  ray  radioactive  sources,  increas¬ 
ing  energy  is  downward;  neutron  residue  increases 
toward  the  right.  The  plots  were  generated  from 
data  taken  from  measurements  of  oscilloscope  traces 
photographed  during  the  experimental  run.  The  numbers 
in  the  plot  refer  to  the  number  of  events  in  each 
"bin".  Source  types  are  indicated  on  each  plot. 

e)  Plot  of  sample  of  raw  neutron  detector  data 
taken  during  run. 

Fig.  13  -  Time  distribution  of  neutron  events  shortly  after 

muon  stopping  time,  taken  with  30  nsec/channel  digi- 
tron.  Early  part  of  data  is  fitted  with  exponentials 
with  decay  times  characteristic  of  Ag,  Zr,  and  Fe 
capture . 

Fig.  14  -  Simplified  block  diagram  of  electronic  logic  used  in 
muon  capture  in  hydrogen  experiment.  The  block 
marked  "logic"  refers  to  the  (i,sec  muon  triggered 
gating  system  used  to  gate  "real"  and  "accidental" 
neutrons  and  electrons. 
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Fig.  15  -  Typical  oscilloscope  traces  photographed  during  run. 
a)  and  b)  Upper  trace  (right  to  left)  shows  a 
neutron  detector  pulse  (negative)  and  a  muon  stopping 
signal  (positive)  and  the  neutron  detector  identifi¬ 
cation  pulse  at  the  left.  Center  trace  shows  the 
neutron  detector  pulse  used  for  energy  measurement 
and  the  slow  residue  pulse  used  for  particles  discrimi¬ 
nation.  The  lower  trace  (high  sweep  speed)  shows  the 
neutron  pulse  and  would  display  pulse  from  detector  #1 
if  it  were  present. 

c)  Upper  trace  shows  an  "accidental"  neutron  event 
due  to  a  second  beam  particle  present  at  the  time  of 
the  neutron  pulse.  The  fast  trace  therefore  displays 
a  pulse  (positive)  from  #1. 

d)  A  negative  pulse  from  the  anticoincidence  detector 
is  visible  on  the  upper  trace. 

Fig.  16  -  a)  and  b)  Time  distribution  of  neutrons  from 

|j,  +  p  -•  n  +  v  •  Time  interval  measurements  were  ob¬ 

tained  from  scanning  of  film;  accidental  events  have 
been  subtracted  but  iron  capture  neutrons  are  in  evi¬ 
dence  at  early  times.  Figs.  16-a  and  16-b  are  inde¬ 
pendent  plots  of  "4  Inch"  and  "3  Inch"  data,  respec¬ 
tively. 

Fig.  17  -  a)  Neutron  pulse  height  distribution;  all  "acci¬ 
dental"  events  from  two  different  neutron  detectors. 
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Fig.  17  -  b)  Time  distribution  of  all  "accidental"  events  as 
measured  on  film.  The  distributions  are  flat  (as 
they  should  be)  except  for  statistical  fluctuations. 

Fig.  18  -  a)  and  b)  Time  distribution  of  neutron  events  occur¬ 
ring  during  positive  beam  running.  Data  from  "4  Inch" 
and  "3  Inch"  geometries  are  presented  separately, 
c)  Pulse  height  distributions  of  neutron  events 
recorded  during  positive  beam  running. 

Fig.  19  -  a)  ,  b)  ,  and  c)  Pulse  height  distribution  of  corrected 
neutron  events  from  p,  +  p  -♦  n  +  v  •  Each  neutron 
detector  and  each  geometrical  arrangement  is  plotted 
individually.  The  horizontal  scale  is  in  MeV  as 
calibrated  with  gamma  ray  sources.  The  results  of  a 
Monte  Carlo  calculation  of  detector  response  are  shown 
as  a  solid  curve.  It  was  fitted  by  normalization. 

Fig.  20  -  Pulse  height  distribution  of  neutrons  from  muon  cap¬ 
ture  in  carbon,  iron,  zirconium  and  silver.  These 
curves  were  measured  in  an  independent  run.  The 
horizontal  scale  is  in  MeV  as  calibrated  with  gamma 
ray  sources.  The  vertical  scale  is  arbitrary  and 
does  not  reflect  the  yield  factors  for  each  of  the 
elements.  The  yields  and  spectrum  shapes  shown  here 
are  not  needed  in  the  analysis  of  the  data  and  were 
only  used  to  estimate  the  effect  of  wall  stopping 


muons . 
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Fig.  21  -  Pulse  height  distribution  of  neutrons  from  muon  cap¬ 
ture  in  iron.  These  points  were  obtained  from  the 
"early  time"  part  of  the  actual  neutron  data  taken 

during  the  run.  The  dashed  curve  shows  how  poor  a 
2 

fit  (x  =  58.3)  is  obtained  if  the  iron  neutrons  are 
fitted  with  the  expected  distribution  from  hydrogen. 

A  much  better  fit  is  obtained  =15.1)  if  an 
exponential  energy  distribution  is  assumed 
[Nq  exp  (-E/0.9) ]. 

Fig.  22  -  The  solid  curve  is  the  pulse  height  spectrum  in  the 
anticoincidence  detectors  due  to  decay  electrons 
which  trigger  an  (NA)  coincidence.  There  are  210 
decay  electrons  contained  in  this  plot.  The  "broken" 
curve  is  the  size  distribution  of  pulses  in  the  anti- 
coincidence  detectors  observed  on  film  during  actual 
data  taking.  The  electronic  anticoincidence  require¬ 
ments  are  in  effect.  In  the  case  of  very  small 
pulses  the  electronic  anticoincidence  fails  and  the 
pulses  appear  on  the  film.  The  pulses  which  fall 
below  the  main  peak  of  the  solid  curve  are  decay 
electrons  which  appear  on  the  film  together  with  an 
"accidental"  neutron  event.  The  "broken"  curve  was 
obtained  from  data  which  corresponds  to  1.2  x  10 
detected  decay  electrons. 

Fig.  23  -  A  simplified  block  diagram  of  the  part  of  the  elec¬ 
tronic  logic  peculiar  to  the  pion  capture  experiment 
(see  Fig.  14  for  notation  and  rest  of  logic).  The 
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pulse  height  analyzer  (PHA)  shown  was  used  to  monitor 

pion  beam  purity.  A  second  PHA  was  gated  on  "events" 

and  displayed  the  neutron  pulse  height  spectrum. 

Fig.  24  -  a)  Gamma  ray  pulse  height  distributions  as  measured 

on  a  PHA  used  to  set  energy  scale  of  neutron  events. 

24  12 

The  sources  are  Na  {2.76  MeV),  C  (4.4  MeV)  and 
016  (6.14  MeV) . 

b)  Uncorrected  pulse  height  distribution  of  neutron 

events  from  pion  capture  (w~  +  p  -*  n  +  7)  .  This 

energy  spectrum  was  taken  on  a  PHA;  it  was  not  used 

in  the  final  analysis  except  as  a  check. 

Fig.  25  -  Pulse  height  distribution  of  pion  capture  neutrons 

taken  from  analysis  of  oscilloscope  photographs. 

The  results  of  a  Monte  Carlo  calculation  of  detector 

response  are  shown  as  a  solid  curve.  The  calculated 

2 

response  is  fitted  by  normalization  and  x  is  com¬ 
puted  afterwards . 

Fig.  26  -  Muon  capture  rates  in  the  (p(i,p)  molecular  ion  with 
the  axial-vector  coupling  constant  and  the  induced 
pseudo-scalar  coupling  constant  as  independent  vari¬ 
ables.  The  point  marked  with  a  capture  rate  of 
560  sec  corresponds  to  the  usual  choice  of  coupling 
constants.  The  band  indicates  the  region  consistent 
with  the  results  of  the  present  experiment. 

Fig.  27  -  a)  The  universal  single  neutron  collision  energy 

2  -A 

response  curve  (A  +  4BE)  showing  spectrum  rise  at 
low  energies  due  to  non-linear  detector  response. 
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Horizontal  scale  is  in  MeV  as  calibrated  with  a 
gamma  ray  source. 

b)  Results  of  a  Monte  Carlo  calculation  for  a  thin 
neutron  detector  with  incident  5  MeV  neutrons  compared 
with  the  theoretical  response  curve.  Actual  response 
measurements  were  also  made  with  a  thin  detector. 

These  fit  the  curve  shown. 

Fig.  28  -  The  pulse  height  response  of  scintillator  NE213  to 

protons.  Experimental  points  were  obtained  using  the 
Nevis  Neutron  Velocity  Selector  (see  text)  and  one 
point  was  obtained  from  pion  capture  measurements . 

The  solid  curves  are  those  of  Brooks  and  Batchelor. 

The  dashed  curve  was  the  one  used  in  the  analysis. 

Fig.  29  -  Radial  distribution  of  muon  trajectories  in  the 

styrofoam  dummy  target  at  several  positions  along 
the  target  axis  (Z) . 

Fig.  30  -  Relative  efficiency  of  the  neutron  detectors  for 
neutrons  originating  within  the  hydrogen  target. 

The  relative  efficiencies  (as  determined  with  a 
Monte  Carlo  calculation)  are  plotted  as  a  function 
of  starting  neutron  energy.  The  several  curves 
show  the  effect  of  the  different  choices  of  a  lower 
energy  threshold  (Ecq) .  The  curves  are  plotted  for 
Ecq  =  0,  1/2,  1,  and  If-  MeV  on  a  gamma  ray  scale. 
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